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Abstract 

In Euclidean relational particle mechanics, only relative times, relative angles and relative separations are mean- 
ingful. Barbour-Bertotti (1982) theory is of this form and can be viewed as a recovery of (a portion of) Newtonian 
mechanics from relational premises. This is of interest in the absolute versus relative motion debate and also shares a 
number of features with the geometrodynamical formulation of general relativity, making it suitable for some modelling 
of the problem of time in quantum gravity. I also study similarity relational particle mechanics ('dynamics of pure 
shape'), in which only relative times, relative angles and ratios of relative separations are meaningful. This I consider 
firstly as it is simpler, particularly in 1 and 2 d, for which the configuration space geometry turns out to be well-known, 
e.g. S 2 for the 'triangleland' (3-particle) case that I consider in detail. Secondly, the similarity model occurs as a 
sub-model within the Euclidean model: that admits a shape-scale split. 

For harmonic oscillator like potentials, similarity triangleland model turns out to have the same mathematics as a 
family of rigid rotor problems, while the Euclidean case turns out to have parallels with the Kepler-Coulomb problem 
in spherical and parabolic coordinates. Previous work on relational mechanics covered cases where the constituent 
subsystems do not exchange relative angular momentum, which is a simplifying (but in some ways undesirable) feature 
paralleling centrality in ordinary mechanics. In this paper I lift this restriction. In each case I reduce the relational 
problem to a standard one, thus obtain various exact, asymptotic and numerical solutions, and then recast these into 
the original mechanical variables for physical interpretation. 
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1 Introduction 



The traditional formulation of mechanics [TJ has rested on absolute space and absolute time. However, Leibniz [2] and 
Mach [3] raised philosophically well-motivated relational objections to these foundations. (See e.g. [H [51 [7J for further 
discussion of this 'absolute versus relative motion debate'.) It is reasonable to consider whether relational principles apply 
to physics as a whole. However, for many years no means were known by which physical theories could be built along 
these lines. Then Barbour and Bertotti [5] and Barbour [9] found some relational particle mechanics [Reissner's earlier 
theory - subsequently rediscovered by Schrodinger and by Barbour-Bertotti - [101 [6] is incompatible with mass-anisotropy 
experiments.] 

Note that the present paper uses the word 'relational' in Barbour's sense. This is worth some discussion because 
Rovelli [TT] uses the same word, but each of he and Barbour take it to mean something different. In outline, Rovelli's 
classical relationalism involves objects not being located in spacetime but being located with respect to each other. Rovelli 
also has a quantum relationalism that quantum states for a particular subsystem only make sense with respect to another 
subsystem. He then speculates that these two relationalisms of his might be related (p 157 of the online version of [TT] : 
"Is there a connection... This is of course very vague, and might lead nowhere, but I find the idea intriguing." Barbour 
on the other hand, has specific spatial and temporal relationalism postulates that embody particular ideas of Mach (that 
time is to be abstracted from change) and Leibniz (the identity of indiscernibles), each of which is sharply implemented 
by particular mathematics at the classical level, as follows. 

A physical theory is temporally relational if there is no meaningful primary notion of time for the whole system thereby 
described (e.g. the universe) [El [321 [13]. This is implemented by using actions that are manifestly reparametrization 
invariant while also being free of extraneous time- related variables [such as Newtonian time or General Relativity (GR)'s 
lapse] . This reparametrization invariance then directly produces primary constraints quadratic in the momenta (such as 
the energy constraint of mechanics or the Hamiltonian constraint of GR). 

A physical theory is configurationally relational if a certain group G of transformations that act on the theory's config- 
uration space Q are physically meaningless [H [12l [Ml [13] . This can be is implemented by such a£] using arbitrary-G- 
frame-corrected quantities rather than 'bare' Q-configurations. For, despite this augmenting Q to the principal bundle 
P(Q, G), variation with respect to each adjoined independent auxiliary G- variable produces a secondary constraint linear 
in the momenta (e.g. the GR momentum constraint arises as a vectorial collection of 3 such constraints) which removes 
one G degree of freedom and one redundant degree of freedom among the Q variables. Thus one ends up dealing with the 
desired reduced configuration space - the quotient space Q/G. Configurational relationalism includes as subcases both 
spatial relationalism and internal relationalism (in the sense of gauge theory). Configurational relationalism can also be 
implemented, at least in some cases [HI [13], by working directly on reduced configuration space. (Reissner's theory was 
formulated along these direct lines, while I also found that the theories [8, 9 can be arrived at in 1- and 2-d from a direct 
implementation [13].) 

One difference between Barbour and Rovelli's approaches is as follows (another is mentioned in Subsec 2.1). In Sec. 
2.4.4 of [TT], Rovelli discusses "meanings of time" and identifies time in Newtonian physics as a metric line that exists 
alongside the configuration space of the system, nowhere reflecting Barbour's starting point that time is derived from 
change alongside consideration of the configuration spaces alone, on which Jacobi-type variational principles are defined. 
This then takes Barbour straight to the situation in which no variable is distinguished as time at the kinematic level; while 
Rovelli characterizes this as an essential difference between non-relativistic and relativistic mechanics in his approach, in 
Barbour's approach this distinction has dissolved. The recovery of Newtonian dynamics from relational particle models 
also has features by which such models are closer in objective structure to GR than Rovelli generally holds nonrelativistic 
mechanics models to be in [TT]. Thereby, relational particle models provide tractable models outside the scheme in [TT] . 

As far as I know, Barbour and collaborators have not as yet reached a specific notion of quantum relationalism. Rovelli's 
own idea of relationalism at the quantum level above-mentioned does share a number of features with record-theoretic 
positions that Barbour and I have previously advocated [T8j [19 , 20] . It is not as yet know whether such records-theoretic 
positions really have any substantial conceptual or technical ties to Barbour's notion of classical relationalism, any more 
than Rovelli knows whether his similar quantum position to be tied to his own classical relationalism (as above-quoted). 

While Barbour's relationalism has the virtue of specifically being in line with Leibniz and Mach (whereby alongside 
having attained a concrete mathematical implementation of these ideas it is definitely of interest to the foundations 
of physics and theoretical physics and so deserves full investigation), I would not dismiss the possibility that Rovelli's 
relationalism is also in line with other (interpretations of) Leibniz, Mach or other such historical figures, and, in any case, 
has original value and is useful in a major quantum gravity program (Loop Quantum Gravity). 

The current paper and its sequel I21| concern technical advances with examples specifically of Barbour's relational 
program - I subsequently use 'relational' in this sense except where I specifically say otherwise. Sec 2 presents Euclidean 
relational particle mechanics [51 [51 [HI [531 (Ml 1251 [Ml EZl H3] (also referred to as scaled models) and similarity relational 
particle mechanics [9, 25, 27] fTTl [13] (also referred to as scalefree models) and further motivates the relational scheme by 
arguing that some (conformo)geometrodynamical formulations of GR can be regarded as arising therein too. Sec 3 explains 
the configuration space structure of relational particle mechanics which has further parallels with GR. All these parallels are 
eventually relevant as regards relational particle mechanics furbishing useful toy models [29l [HI [191 ESI [271 EE EQ1 El 132] 

lr This is my own 14 15 16 passive, mathematician's implementation, whereas Barbour thinks about this in active terms that physicists 
sometimes use. This difference in thinking does not, however, lead to any tangible discrepancies in the material of this paper. 
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for the Problem of Time in Quantum Gravity [551 133] and other issues in the foundations of Quantum Cosmology [TH1 [31] • 
Use of relational particle mechanics in both the absolute versus relative motion debate and the study of conceptual 
strategies suggested toward resolving the Problem of Time in Quantum Gravity would benefit from having a good working 
understanding of explicit examples of relational particle mechanics, which is the subject of this paper at the classical level. 
Relational particle models concern N particles in dimension d. As the general configuration for d = 3 is an N-haedron, 

1 term the 3-d N-particle relational particle model N-haedronland. Likewise, I term the 2-d N-particle relational particle 
model N-a-gonland, and the 1-d one N-stop metroland (as in urban public transport maps). The special 3-a-gonland 
considered as the principal example in this paper I term triangleland. 

Noting that scalefree N-stop metroland and N-a-gonland have fairly standard geometry (§ N_2 and CP N_2 respectively 
|13| ) permits explicit reductions [171 113] and subsequent availability of useful coordinate systems and methods of math- 
ematical physics. For scalefree triangleland, CP 1 = § 2 , so one has 'twice as many techniques', so in this paper I choose 
to study this case. Scaled triangleland also permits explicit reduction and its configuration space in shape-scale variables 
takes the form C(§ 2 ), which also makes for tractable and interesting explicit examples, where the cone [35] C(X) over a 
space X is R+ x X [J 0, the special cone-point or 'apex'. 

In Sec 4 I give Euler-Lagrange equations for scalefree N-stop metroland and N-a-gonland and give further specific 
forms for the exceptional triangleland case. I then consider examples with potentials that are independent of the relative 
angle <l|l as well as harder ones that depend on this. (The former is a substantial simplification [T7] - in close analogy 
with centrality in ordinary mechanics - but it is the harder latter case that is relevant to various Problem of Time schemes 
- semiclassical emergent time and possibly the semblance of dynamics in timeless records schemes discussed in paper II 
[21|.) In particular, I consider the general case of harmonic oscillator like potentials between all particles, looking at the 
'special' (^-independent) and 'very special' (constant) subcases within as well as the general small and large asymptotic 
behaviour, and recast this in terms of the problem's 'original variables' - mass-weighted relative Jacobi variables i x , i 2 
that still contain an absolute orientation. Identifying the 'very special' problem's mathematics as corresponding to the 
linear rigid rotor, the 'special' problem's as that with additionally a background homogeneous electric field in the axial 
direction and the general problem's likewise but now with a general direction, is valuable in this paper and in subsequent 
work at the QM level in paper II. In Sec 5 I use that a rotation (or, cquivalently, a normal modes construction) maps the 
general case to the special case for this particular problem. Thus I can get as far in solving for the general case as I can 
with the 'special' case. 

In Sec 6, 1 consider scaled triangleland's Euler-Lagrange equations in terms of the straightforwardly relational variables 
(ti,i 2 )<&), the useful (Ii,I 2 ,$) coordinates (that turn out to be parabolic coordinates), the C^CP 1 ) presentation shape- 
scale coordinates {I,TZ, <&) and the C(§ 2 ) presentation shape-scale coordinates (I, 9, $). I then exactly solve the 'special' 
case of multi-harmonic oscillator like potential for scaled triangleland, by mapping it in (Ii, 72,3?) coordinates to a close 
analogue of the Kepler-Coulomb problem (which move remains useful at the quantum level 37J). I then give Euclidean 
relational particle mechanics' own rotation/normal modes construction, whereby my having obtained the special solution 
enables me to also obtain the general solution, albeit this paper only has room for presenting the special solution and its 
physical interpretation. 

I conclude in Sec 7, including an outline of further promising relational particle mechanics examples. Paper II considers 
this paper's similarity models at the quantum level and 37J likewise for this paper's Euclidean models. Interesting Problem 
of Time in Quantum Gravity applications of these models will be developed in yet further papers [31] . 

2 Examples of relational theories 
2.1 Euclidean relational particle mechanics 

In Euclidean, or scaled, relational particle mechanics [51 [51 [2"31 119] . only relative times, relative angles and relative sep- 
arations are meaningful. E.g. for 3 particles in dimension d > 1, Euclidean relational particle mechanics is a dynam- 
ics of the triangle that the 3 particles form. Euclidean relational particle mechanics was originally [8J conceived for 
Q = Q(N, d) = R Nd the positions q Ia of N particles in d-dimcnsional space with G the d-dimensional Euclidean group 
Eucl(d) of translations, Tr(d), and rotations, Rot(d)H However, eliminating Tr(d) is trivial and produces a theory of 
essentially the same form as the original if relative Jacobi coordinates Ri a are employed |26j , so I take that as my starting 
position. Relative Jacobi coordinates are inter-particle (cluster) separations chosen such that the kinetic term is diagonal 

2 The notation for this and the next paragraph is as follows. R^, i = 1, 2, are relative Jacobi coordinates [36]. 3> is the angle between 
these. By being mass- weighted, I mean that t i = ^/JIiR it where fij are the particle (cluster) masses associated with R i 1361 . (.j = 
the magnitude of i i , and 1% = l\ = p^iR\ , the ith Jacobi (barycentric) partial moment of inertia. 1Z is the simple ratio variable 11/12 and 
© = 2arctan7£, which turns out to geometrically be the azimuthal spherical angle; see also Fig 1 and Sees 2—3 for further interpretation and 
depictions of these. 

3 Lower-case Greek letters are spatial indices. Capital indices are used for the N particle position coordinates, lower-case indices for n = N — 
1 coordinates describing relative particle (cluster) separations, barred and tilded lower-case indices to take values 1 to n - 1 values and hatted 
lower-case indices to take values 1 to n - 2. The dot denotes d/dA for A a label-time parameter that has no physical meaning since JTJ is 
reparametrization-invariant. B a generates the rotations. As well as the obvious 3-d case, the 1 and 2 d cases are incorporable into this form as 
follows. Take B a = (0, 0, B) in 2-d. Then L a has just one component that is nontrivially zero. Take furthermore B = in 1-d. Then there is 
no L" constraint at all. M IJa ^ is the kinetic metric )j,i5 l: >5 a P , where fn are the Jacobi masses |36| . [E.g. for this paper's triangleland example, 
these are fi± = m^ra-j, / '{m.2 + m^} and fj.2 = m\{m,2 + m:j}/{mi + + "13} for mj the particle masses.]. J^ij a p is the inverse of this array. 
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[55] : thus I take my Q to be R(N, d), the space of relative separations (this is R nd ), and G to be Rot(d). Then configu- 
rational relationalism takes the form that one is to construct one's action using the arbitrary Rot(d) frame expressions 
oR ia = R ia — tcP" 1 BpRi 1 rather than 'bare' Ri a . 

[Figure 1: Coordinate systems for scaled triangleland. 

i) Absolute particle position coordinates [q_^, <j 2 , q^) with respect to fixed axes and a fixed origin O; the corresponding particle masses are mj, 
I = 1 to 3. 

ii) Relative particle position coordinates, any 2 of which form a basis. 

iii) Relative Jacobi coordinates (R_i, R?)', the corresponding Jacobi masses are fa, i = 1, 2. The mass-scaled relative Jacobi coordinates are 
related to these by l ( = ^/JmBh- 

iv) Bipolar relative Jacobi coordinates (pi, 9\, pi, Bi). The mass-scaled radial Jacobi coordinates are n = sjjipi- These coordinates still refer 
to fixed axes. 

(R -R \ 
iTg^Trfjjl || J th e relational 'Swiss army knife angle' between the 2 relative Jacobi 

vectors. The coordinate ranges are 0<pi<oo,0<$< 2n. ] 

The action that one builds is of Jacobi type |38j so as to implement temporal relationalism, 

S[R ia ,R ia ,B a }=2 J dAv/T{U + E} . (1) 

Here the kinetic term T(R ia , R ia , B a ) = M l i a P o R ia o Rjp/2 is homogeneous quadratic in the velocities. U is minus the 
potential energy V which is a function of y/Rj, ■ Rj alone, and E is the total energy of the closed system/universe; hitherto 
this has been taken to have a fixed value (but should it: no observer would know it exactly?) [Note that each such action 
Sj = 2 J dA^/TjE + U} is indeed equivalent to the more well known Euler-Lagrange actions 



S = J dt{T - V} , (2) 



where T now differs from the its previous form in containing * in place of '. See e.g. |38j for obtaining actions of type ((T|) 
from actions of type ([2]) by parametrization and Routhian reduction to eliminate di/dA.] 
From ([I]), the conjugate momenta are then 



P ia = M ljal3 * Rj$ for * R la = Riot - e a ^Bf 3 *R ll (3) 



and * = a/{U + E}/T ' = d/dt for t the emergent 'Leibniz-Mach-Barbour' time that coincides here with Newtonian time. 
This object introduced, another difference between Barbour's scheme and Rovelli's can be pointed out. In Rovelli's scheme, 
one clock/timestandard is much as good as another (at the conceptual level rather than what is convenient or accurate), 
while in Barbour's scheme there emerges a Leibniz-Mach-Barbour timestandard to which everything in the universe 
contributes (an 'ephemeris' timestandard [S3J 133 SO] ) ■ This choice of time is distinguished by its substantially simplifying 
both the above momentum- velocity relations and the Euler-Lagrange equations, and it amounts to an emergent recovery 
of other notions of time like Newtonian time here in mechanics, or proper time or cosmic time in Sec 2.3, and it has some 
parallels with the actual ephemeris time in official use in the first half of the 20th century [39] ■ which is based on the 
totality of the motions of the objects in the solar system. 

Reparametrization invariance implies [41j that these must obey at least one primary constraint; here there is one, 
which takes the form of an energy constraint 

H = iV yQ/3 P M P^/2 + V= E , (4) 

to which the momenta contribute quadratically but not linearly. Variation with respect to B a yields as a secondary 
constraint the zero total angular momentum constraint 

L" = e a ^R lfj P^ = , (5) 

which is linear in the momenta and is interpretable as the physical content of the theory being in relative separations and 
relative angles and not in absolute angles. 

Euclidean relational particle mechanics is Leibnizian/Machian in form, and yet is in agreement with a subset of 
Newtonian mechanics - the zero total angular momentum universes. The recovery from relationalism of many of the results 
previously obtained by standard absolutist physics and the finding of similar mathematical structures along both routes 
(at least in the simpler cases) is an interesting reconciliation as regards the extent to which we have not been prejudiced by 
studying absolutist physics. This paper and [5T] show that this trend extends further, albeit the mathematics arising from 
similarity relational particle mechanics is distinct from that in the ordinary central force problem in having restricted, 
unusual potentials inherited from the scale invariance). Sec 2.3 and 2.4 recollect that a similar trend is also present in 
GR. 



3 



2.2 Similarity relational particle mechanics 



In similarity, or scalefree, relational particle mechanics [51 1271 [T7] . only relative times, relative angles and ratios of relative 
separations are meaningful. I.e., it is a dynamics of shape excluding size: a dynamics of pure shape. In the case of 3 
particles in dimension d > 1, similarity relational particle mechanics is the dynamics of the shape of the triangle that the 
3 particles form. G is now augmented by the dilations Dil(d) = Dil to form the Similarity group Sim(d), so that oR ia 
now takes the form oR ia = R ia — e a ^ 'BpR^ + CRi a , where C generates the dilations. Noting that the 'banal conformal 
transformation'0 

T — * Tn = ^ 2 T , E + U^E n + U f2 = {E-V}/f2 2 , (6) 
leaves the Jacobi action invariant, the most natural presentation [17) for the action is 

S[R ia ,R ia ,B a ,C]=2 J dAv/T{U + E} (7) 

with T(Ri a , Ri a , B a , C) = M u "^ o R ia oRj p/21 for I the total barycentric moment of inertia of the system, whereupon V is 
a function only of manifestly scale-invariant ratios of yjR[~Rj (i.e. homogeneous of degree zero) and E comes unweighted. 
[Barbour's original presentation [S], in which the potential is homogeneous of degree —2 and an energy E cannot be added 
on to this (but E/I can), is related to this one by use of SI 2 = I, i.e. by not dividing through by the moment of inertia 
in the first place.] Then (in the conformally natural presentation I use above) the conjugate momenta are given by ([3]) 
but divided by I and containing the * corresponding to the new o. There are again primary and secondary constraints 
respectively of form ^ with the first term multiplied by I and ([H]), as well as a secondary constraint from variation with 
respect to C, 

D = R la P la = , (8) 

which is also linear in the momenta. This is the dilational (or Euler) constraint - it says that the dilational momentum of 
the whole system is zero so that the physical content of the theory is not in relative separations but in ratios of relative 
separations (relative angles already being functions of ratios, they are unaffected) . 



2.3 Relational formulation of geometrodynamics 

Important further motivation for the study of relationalism and relational particle mechanics is that General Relativity 
(GR) can also be formulated as a relational theory. Begin by recollecting that as well as being a spacetime theory, GR can 
be studied as a dynamics obtained by splitting spacetime with respect to a family of spatial hypersurfaces [45] . However, 
answering a question of Wheeler [46] , this dynamics can be taken to follow from first principles of its own [47l [12l Q~5] . 
And one of the two known such sets of first principles are relational first principles, with Q = Riem(E, 3) - the space of 
positive-definite 3-metrics on some fixed topology £ which I take to be a compact without boundary one for simplicity - 
and G = Diff(E, 3) the diffeomorphisms on £ (TSJ [14l [15] - The arbitrary Diff(S, 3) frame expressions are 
o/i M „ = h^ u — £ ph^ rather than 'bare' /i^H 

The action that one builds so as to explicitly implement temporal relationalism is [T21 HH HPJ US] 

S GK [h ap , h af3 ,F ] = 2 J AX J d 3 xVhyjT GK {Ric(h)-2A} for T GR (z*\ V> V; K\ = \^ vp ° o V o h pa , (9) 

which bears many similarities to the better-known Baierlein-Sharp- Wheeler [J5] action@ Also note that one does not 
need to assume the GR form of the kinetic metric or potential; relational postulates plus a few simplicities give this since 
the Dirac procedure [H] prevents most other likewise simple choices of kinetic term T from working [I^J [SOI EH Q31 Q~5] . 
Yet further motivation is that 1) configurational relationalism is closely related [T3JI23] to certain formulations of gauge 
theory. 2) The above relational formulation of GR is furthermore robust to the inclusion of a sufficiently full set of 
fundamental matter sources so as to describe nature [52[ [53l [14l [T5j [16] . 

4 Performing this transformation clearly leaves invariant product- type actions and so of Jacobi- type actions JTJ and its similarity relational 
particle mechanics counterpart and of the reduced actions that follow from these in Sec 3, and of their GR counterparts such as the Baierlein- 
Sharp- Wheeler action 42 , J9j, and J 13D and its variants. T and E + U scale compensatingly and then deduce from * = ^/{E + U}/T ' that * 

scales as * > *Si = n~ 2 *. (Euler-La grange or Arnowitt-Deser— Misner type actions have this invariance too but its presence therein is less 

obvious to spot 43 .) I term each particular choice of SI a 'banal conformal representation'. Note: classically Euler-Lagrange equations are 
invariant |43j so applying a banal conformal transformation makes no difference and it is just a means of computational convenience in some 
cases. But retaining this lack of difference at the quantum level has implications |43II21||3"T1 . 

5 The dot now denotes d/dX. F^ generates Diff(S, 3). £p is the Lie derivative with respect to the vector field _F M . h, and Ric(/i) are 
the determinant, covariant derivative and Ricci scalar associated with A is the cosmological constant. A4'" /p<T = Wh 1 " 7 — h^h pa is the 

kinetic supermetric of GR, with determinant M and inverse N^vop = hn P h va — h^hpu /2 (which is the undensitized version of the DeWitt 
supcrmetric [48]). 

6 The difference is that the Baierlein-Sharp-Wheeler action contains the shift, (3^, while the action J9} contains the velocity associated with 
the frame variable Fu,, which is such that F M = Both of these actions are free of extraneous time-related variables (unlike the Arnowitt- 
Deser-Misner action which contains a such — the lapse - multiplier elimination of which produces the Baierlein-Sharp-Wheeler action). The 
above-mentioned difference does not affect the outcome of the variational procedure |44| ; however, it does make |(9j homogeneous quadratic in 
d/dX so that the A's cancel out of dA\/T so that it and not the Baierlein-Sharp— Wheeler action is manifestly reparametrization invariant as 
well as free of extraneous time-related variables, i.e. temporally relational as defined on page 1. 
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Then the conjugate momenta are 



tt^ = sfhM^ vp<x * h pa for * h pa = V* - £ F *h pa (10) 

for * = \J {Ric(ft-) — 2A}/T GR ' . ([9]) being reparametrization invariant, there must likewise be at least one primary 
constraint, which is in this case the GR Hamiltonian constraint 

H = -^A/^p (T 7T^7r p<T - Vh{Ric{h) - 2A} = (11) 
Vh 

to which the momenta contribute quadratically but not linearly. Variation with respect to F p yields as a secondary 
constraint the GR momentum constraint 

H p = -2A,tt% = , (12) 

which is linear in the gravitational momenta and interpretablc as GR being more than just a theory of dynamical 3-metrics: 
the physical information is in the underlying geometry and not in the allocation of points to that geometry. While, the 
purely quadratic nature of Ti leads to ji^l = and not ihd^/dT for some notion of time T. The H of relational particle 
mechanics also has this feature and so itself manifests the Problem of Time (Sec II. 1); it is a reasonable model for this 
in that a number of the conceptual strategies subsequently suggested for GR have nontrivial counterparts for relational 
particle mechanics. 

The presence of square roots in the above actions so as to implement manifest temporal relationalism has caused some 
concern with Referees due to square roots causing substantial difficulties at the quantum level. In this respect, I note that 
the above square roots occur at the classical level in the Lagrangian formulation. However, there are no square roots in 
the Hamiltonians resulting from such actions, and it is these that I promote to quantum equations in 21j, so that there 
are no square roots in the quantum equations, and so this work encounters none of the problems associated with handling 
square roots at the quantum level. 

2.4 Relational formulation of conformogeometrodynamics 

This has a number of ties with the work of Lichnerowicz [54j and York [55] on maximal and constant mean curvature 
spatial slices respectively, which is important in numerical relativity [56] . A preliminary action is 

S GR = J dxj d 3 xVh<t> 6 ^T GR {(t>- 4 {Ric(h) - 8D 2 (f>/^}} , (13) 



wher^l 



{WW - h^h^}{oh^ + Ah pv o ^/4>}{oh pa + Ah p<J o (/,/</,} ■ (14) 
is a conformal factor. This action then gives as a primary constraint 

12, 



= t-s-Lw^ -Vh(mc{h) - = o 



(15) 



(Lichnerowicz equation), TL p as a secondary constraint from F p variation, and 

7T = (16) 

(condition for a maximal slice) as one part of the free end hypersurface [4"4"l [55] variation of <f>, but the other part of this 
variation entails frozenness (the well-known non-propagability of maximal slicing for spatially compact without boundary 
GR. 

[57] and |49j got round this frozenness was circumvented by considering a new action with division by Vol 2 / 3 which 
amounts to fully using G = Diff(S, 3) x Conf(E, 3) where Conf(S, 3) is the group of conformal transformations on S. 
But the subsequent variational principle no longer gives GR and is furthermore questionable as an alternative theory 
[491 El [51]). 

In [58] it was circumvented rather by using G = Diff(E, 3) x VPConf(E, 3), where VPConf(E, 3) are the (global) 
volume-preserving conformal transformations on S as implemented by using 

1 1/6 

d 3 xVh(j) 6 \ (17) 



'I, like |57) . present this for simplicity with no A term; see [49] for inclusion of this. Also, |57| first considered this with two separate 
multipliers instead of a single more general auxiliary whose velocity also features in the action and has to be free end hypersurface varied; the 
way this is presented here is that of |58j and 1591 . See also |44l 1591 for justification of the type of variation in use. 
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Subsequently, the primary constraint (now denoted by H?) picks up more terms (making it a Lichnerowicz-York equation 
rather than a Lichnerowicz equation). While, 4> variation now gives 

ir/Vh = C, (18) 

(condition for a constant mean curvature slice) alongside an equation that successfully maintains this. The addition of 
matter to the present paragraph's scheme has not to date been extensively studied, but no significant hindrances are 
known to date [SD] , and the addition of matter to the preceding paragraph's scheme has been extensively studied [121 [M] • 



3 Geometry of the configuration spaces 

The following configuration space structure issues play an important underlying role in papers I and II. 



3.1 Euclidean relational particle mechanics configuration spaces and reduction 

For Euclidean relational particle mechanics, such a study was carried out in [24| in a 'rigged' fashion (auxiliary variables 
included), and in |621 1241 126j in a reduced fashion (auxiliary variables eliminated). [63] exemplifies related work. In 
[2"7] I noted that this elimination has a simpler nature in 2-d than in 3-d, and used this to pass to completely relational 
variables in |17H13j for both Euclidean and similarity relational particle mechanics. Relative space R(N,d) is the quotient 
space Q(N, d)/Tr(d). This is flat space, so the kinetic metric is just T(Ri) = J2 i HiRf /2. At this level have Jacobi 
coordinates and constraints as in Sec 2.1. I use mass-weighted Jacobi variables as the most succinct 'original variables of 
the problem'. Relational space is 7£(N, d) = Q(N, d)/Eucl(d), which is what one has on eliminating the angular momentum 
constraint. These coincide in 1-d. In 2-d one can do elimination explicitly [T71II3], obtaining e.g. in the 3-particle case 
the T tnange (ij, Li, $) corresponding to the line element 

d4iangle = ^\ + + ^T~^ ' ( 19 ) 

h. L 2 

[To obtain a configuration space kinetic term T from the line element of a metric ds 2 = MABdQ A dQ B , use (|155|) .] 
This elimination is done by: Step 1) passing to mass-weighted Jacobi coordinates. Step 2: eliminating the rotational 
auxiliary from the Lagrangian form of the constraint l|5|). casting the subsequent expression in mass- weighted Jacobi 
bipolar coordinates, which I denote by (n, 9\, 12, 62) [fig liv); these are still with respect to fixed axes.] Step 3: one 
can then pass to fully Euclideanly-relational coordinates [fig lv)] ({ti,i2,$}) (by an additional 'purely absolute' angle 
dropping out of the working [26[ I17j ). In 3-d, the situation is considerably harder. 



3.2 Similarity relational particle mechanics configuration spaces and reduction 

For similarity relational particle mechanics, such a study was carried out in 13] in a reduced fashion (both by reduction 
and in an already-reduced form based on [64]). Preshape space is P(N,d) = Q(N, d)/Tr(d) x Dil and Shape space is 
S(N, d) = Q(N, d)/Sim(d). Furthermore, P(N, d)=S nd_1 and S(N, l)=§ n - x where in both cases '=' here means equal both 
topologically and metrically, with the standard (hyper)spherical metric. Thus one has T N_stop (7?.p, IZp) or T N_stop (9p, Op) 
corresponding to the line element 



{1 + n} - {e^i 1 n p <m p ) 



n— 1 r—1 



d 4-sto P = i J n , ^n-i^ 2 ; 2 = E II sin 2 epd9 2 (20) 



r=l p=l 



in terms of simple ratio coordinates (|21|) and (ultra) spherical coordinates (|22p. respectively, where Yli=i t erm s are defined 
to be 1. It is obtained from the similarity relational particle mechanics action in Jacobi coordinates by: Step 1) passing 
to mass-weighted coordinates. Step 2) eliminate dilational auxiliaries from the Lagrangian form of the constraint © and 
re-express in mass-weighted Jacobi bipolar coordinates. Step 3) pass to fully relational simple ratio variables 

Tlp = ip/in-l (21) 

(in the present context these are, geometrically, Beltrami coordinates); these are related to ultraspherical coordinates by 



Qp — arctan 



j>^ P - + i] (22) 



for 72-n-i = tn— i/tn — 1 = 1. The coordinate ranges for these are Qp G (0, 7r), 9 n d-i € [0, 2n). E.g. for 4-stop metroland 
0i is the azimuthal angle and 2 is the polar angle $. Then one has T 4 - st °P(7^ 7^.) or T N - Stop (9,0,$) corresponding 
to the line element 

d.u P - 1 1 + Ri!dRf + g ; g^g F _ ^ + sinW (2 3 ) 



G 



for 'simple ratio' coordinates (IZi, IZ2) [4-stop metroland subcase of (|21[) ] and spherical coordinates (9, $) [ = (Oi, 82), 
given by the 4-stop metroland subcase of (|22|) ]. 

While, for N particles in 2-d S(N, 2)=CP n_1 both topologically and metrically, with the standard Fubini-Study metric 
corresponding to the line element 

, 2 = {1 + Eg \z P \ 2 } E, \Az q \ 2 - I S P -^dZ p -| 2 

^N-a-gon {1+Efl^r| 2 } 2 ' 

for Z F inhomogeneous coordinates on CP 11 " 1 , from which the kinetic term T N_a ~ s ° n (.2p, Zp) is constructed. This is 
obtained from the similarity relational particle mechanics action in Jacobi coordinates by Step 1) passing to mass- weighted 
coordinates. Step 2) eliminate both a rotational auxiliary and a dilational auxiliary from the Lagrangian forms of (0 [8J) - 
Step 3) write this action in terms of ratio variables. Furthermore, one can use the polar form Z = 7Zpexp(i<dp); indexing 
moduli (real ratio coordinates) as IZp and arguments (relative angle coordinates) as Op (for both p and p taking 1 to N - 
2). Then the configuration space metric can be written in two blocks {Mpq = 0): 

Mp 9 = {l + \\K\\ 2 }- 1 8 i %-{l + \\K\\ 2 }- 2 1%Tl 9 , M M = {{l + H^I^-^^-ll + ll^ll 2 }- 2 ^^}^^ (no sum) , (25) 

where for a given p IZp is the IZp that forms a complex coordinate pair with Qp. 

As a particular example, S(3, 2) = CP 1 = § 2 , making this triangleland case particularly amenable to study due 
to the availability of both projective and spherical techniques (while, as we shall see, this example meets many of the 
nontrivialities required by Problem of Time strategies); this paper principally considers this case. In this case the kinetic 
term collapses to T tllan s lc (j7, Z), T^^ 11 ^ (1Z, 1Z, 6) or T|" 2 a "^(8, 9, $) as constructed from the line element 

ds 2 ri a„ g ic = |dZ| 2 /{l + |Z| 2 } 2 = {&1Z 2 + ft 2 d$ 2 }/{l + ft 2 } 2 = {d9 2 + sin 2 9d$ 2 }/4 , (26) 

where TZ is the simple ratio variable TZ = (4/12 This is physically the square root of the ratio of partial moments of inertia, 
■y/li/12, and mathematically a choice of inhomogeneous coordinate's modulus on CP 1 . This is related to the azimuthal 
angle 9 on § 2 = CP 1 by 

9 = 2arctanft K = tanf (27) 

for 9 the spherical azimuthal angle (see also Fig 2). [Thus 1Z is also geometrically interpretable as a radial stereographic 
coordinate on § 2 ]. $ is the relative angle between the two Jacobi coordinate vectors. Useful relations between the two 
representations include that 1Z = 1 is the equator, 1Z small (compared to 1) corresponds to near the North Pole (9 a small 
angle) and 1Z large corresponds to near the South Pole (the supplement angle 3 = it — 9 is small). Also, one can use the 
barred banal conformal representation T — > 4T, U + E — > 4{U + E}, whereupon this sphere of radius 1/2 becomes the 

unit sphere, and I write T§" anslc . 

[Figure 2: Interrelation between TZ, U, © and S coordinates. N is the North Pole, S is the South Pole, O is the centre. Then point P in the 
spherical representation corresponds to point P 1 in the stereographic tangent plane at N with radial coordinate 7?., and to point P" in the 
stereographic tangent plane at S with radial coordinate U.\ 



For use in paper II, this line element in terms of $ and F, i = 1 or 2 is (no sum) 

As 2 = dI 2 /4F{I - F} + {I - F}Fd$ 2 /I 2 . (28) 

— triangle • 

One can also use Tfl a t = T (1Z,1Z, 3>) constructed from the line clement 

d s 2 at =7^ 2 +ft 2 $ 2 (29) 

by performing a banal conformal transformation with conformal factor f2 2 = {1 + 7Z 2 } 2 . This is geometrically trivial, 
while the other above forms are both geometrically natural and mechanically natural (equivalent to E appearing as an 
eigenvalue free of weight function); using {1 + 7Z 2 } alone would be the conformally-natural choice. 

In 3-d, the situation is, again, harder 64J though at least similarity relational particle mechanics is free of collisions 
that are maximal (i.e. between all the particles at once). 



3.3 Euclidean relational particle mechanics in scale— shape variables 

Finally, these 'shapes' are also relevant within the corresponding Euclidean relational particle mechanics, as these admit 
conceptually interesting formulations in terms of scale-shape split variables. For scaled N-stop metroland, the configu- 
ration space is the generalized cone CiS 11 - 1 ) while for N-a-gonland it is CXCP 1 ^ 1 ), scaled triangleland also being C(§ 2 ). 
The special cone-point or 'apex' here corresponds physically to the maximal collision. Now use the new coordinate 

I = l\ + i\ , (30) 
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which is physically the moment of inertia and mathematically a radius, and the same TZ coordinate as in Sec 3.2. This 
coordinate transformation then inverts to 

41 = VW+^jK , 4 2 = y/l/{l+K 2 } . (31) 

Also, TZ can be supplanted by 8 = 2arctan7?.. 

The kinetic term is T trianglo (I, Z, I, Z), T trianglc (I, TZ, I, TZ, $) or T trianglc (I, 9, 1, 9, $) as constructed from the line 
element 

d4 iangl e ^ | {dl 2 + 4I 2 TT ^_} - £ {dl 2 + = k ^ + l2{dQ2 + Sin2ed<&2 ^ • (32) 

Moreover, one can use instead the banal-conformally related kinetic term TFi at (I, 9, 1, 9, $) constructed from 

ds| lat = dl 2 + I 2 {d9 2 + sin 2 9d$ 2 } (33) 

(the corresponding conformal factor being 41, so that also 

U + E = {U + E}/4l) , (34) 

away from I = in which place this conformal transformation is invalid). The form of (|33[) clearly makes this a flat (i.e. 
geometrically trivial) representation in spherical polar coordinates, with I as the radius. While, the other forms above are 
mechanically natural]! 



3.4 Geometrodynamical configuration spaces 

While, GR is a geometrodynamics [351 [33] on the quotient configuration space Superspace(S, 3) = 

Riem(£, 3)/Diff(£, 3) [35J[35], which is studied topologically and geometrically in e.g. [35J [35J [55] . Superspace is an 
infinite-dimensional complicatedly stratified manifold; explicit reduction is not in general possible here. 



3.5 Conformogeometrodynamical configuration spaces 

One can obtain relational theories or formulations using G = Diff(E, 3) x Conf(S, 3) that reproduces the maximal 
condition (however this formulation freezes unless one alters one's theory from GR [49]) and using G = Diff(£, 3) x 
VPConf(E, 3) (the volume-preserving conformal transformations) that does reproduce GR in the York formulation from 
an action principle 58J. 

Here the associated configuration spaces are conformal superspace CS(E, 3) = Riem(£, 3) /Diff(E, 3) x Conf(£, 3), 
which is studied geometrically in [55J[57], and {CS + V}(£, 3) = Riem(£, 3)/Diff(E, 3) x VPConf(E, 3), which previously 
featured in e.g. [55J, but has not to my knowledge been studied from a geometrical perspective. [The CRiem(E, 3) = 
Riem(E, 3)/Conf(E, 3) analogue of preshape space has been studied geometrically in e.g. [48] [67] ]. 

Parallels between this and Euclidean relational particle mechanics in shape-scale split variables include homogeneous 
degrees of freedom playing a similar role to scale, York time [55] having an 'Euler time' analogue [57], similarity relational 
particle mechanics corresponding to maximal slicing in having this time frozen, and in the analogy between {CS + V}(£, 
3) and the Euclidean relational particle mechanics configuration spaces' cone structure. In the configuration space of GR, 
the analogous 'special point' at zero scale is the Big Bang. 



4 Similarity relational particle mechanics at the classical level 
4.1 1 and 2-d cases 

In App A, I provide how the general curved space mechanics unfolds. In the case of scalefree N-stop metroland, the Jacobi 
action is 

S N - stop (9 p -, 9 P -) = 2 j dAv/Tsn-^E-V} . (35) 

The Euler-Lagrange equations following from this are given in [13j in Beltrami coordinates (the notation there uses Sp 
in place of the present paper's TZp), but in this paper I use, rather, (ultra)spherical coordinates in the present paper, so I 
provide the Euler Lagrange equations afresh in these: 



'q-l 1 f n-l f-1 

rising - e IT 

P J I r=q+l p=l^q 



\{ sin 2 9p9;- \ - \ II sin2Q P \ sin9 g -cos9 9 -9; 2 = . (36) 



s 1351 connects this section's workings with results from the celestial mechanics and molecular physics literatures. 
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There is also a first energy integral, 

n- 1 f-i 

II sin2 ep6* 2 /2 + V(0 F ) = E . (37) 

r—1 p—1 

For scalefree N-a-gonland, the Jacobi action 

S = 2 J dAv/T FS {E-V} (38) 

gives the Euler-Lagrange equations as presented implicitly in |13j . 

In the exceptional triangleland case, there are spherical flat (conformal-to-stereographic) presentations, 



5 =2 J dA^/T s2 {E - V} = 2 J dA^WE - V} . (39) 
Then, in the spherical presentation, the Euler-Lagrange equations simplify to 

9** - sinecos$$* 2 = -- , {sin 2 6$*}* = (40) 
with an accompanying energy integral 

9 T2 /2 + sin 2 6$* 2 /2 + V = E . (41) 

[Above, 



'{E + U}/T • = V{E + U}/T/4 • = */4 .] (42) 
If V is independent of $, then (f40l i) becomes another first integral: 

sin 2 9$ 7 = J , constant . (43) 
This is a relative angular momentum quantity; see App C for various interpretations of it. Then (|4ip becomes 

9*2/2 + t 7 2 /2sin 2 e + V(6) = E . (44) 
Or, in terms of the TZ coordinate, the important results (|4"3"|) and (|4~4"1) take the form 

^= j -\\{ish J A= i+ "' (45) 



where * = y {E + U}/T'. This is illuminating though its direct parallel with the usual flat planar presentation of the 
ordinary mechanics of a test particle moving in a central potential: 

ratio of square roots of the two subsystems' \ 

Jacobi partial moments of inertia J = TZ < — ► r = (radial coordinate of test particle) , (46) 

(relative angle between Jacobi coordinates) = $ < — ► 8 = (polar coordinate of test particle) , (47) 

(relative angular momentumX / angular momentum componentX 

of the two subsystems = J < — ► L ( = L z ) = perpendicular to the plane , (48) 

1 < — ► m = (test particle mass) . (49) 

This analogy is furthermore a pointer to parallel the well-known u = 1/r substitution by the hi = 1/TZ one, which turns 
out to be exceedingly useful in my study below. [In the spherical presentation the counterpart of this inversion map 
TZ — ► hi = 1/TZ is the supplementary map — ► tt — = the supplementary angle, S. In or Ii variables, it takes the 
form of interchanging the 1-indices and the 2-indices. I term the underlying operation that takes these forms in these 
presentations as the duality map.] 

Next, V e ff = V + J 2 /TZ 2 — E is the potential quantity that is significant for motion in time, and combining (14"4"1) and 
(|43"]) . U or b = — 7?. 4 V ff is the potential quantity that is significant as regards the shapes of the classical orbits. Translating 
into the spherical language, the corresponding quantities are V c ff = V + J 2 /sin 2 — E and U rb = — sin 4 0V o ff. 

Finally, combining (j44| and (|43|) or their 7?.-analogues gives as quadratures for the shapes of the orbits 

$ - $o = J JdQ/smQ^I 2{E - V(6)}sin 2 e - J 2 = J JdTZ/TZ^2{E + U}TZ 2 - J 2 . (50) 
While, ([33]) and its TZ analogue give as quadratures for the time-traversals 

t-t = J sinede/y^E + Ulsn^e-./ 2 , or t-to = J TZATZ/\j2{E + (i}TZ 2 - J 2 . (51) 
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4.2 A class of separable potentials for scalefree triangleland 

As 9- independent (central) potentials considerably simplify classical and quantum mechanics, by the above analogy it 
is clear that ^-independent (relative angle independent) potentials will also considerably simplify classical and quantum 
scalefree triangleland. In particular, these simplifications include separability. Also, in each case, there is a conserved 
quantity: angular momentum L in the case of ordinary mechanics with central potentials, and relative angular momentum 
J of the two constituent subsystems in the present context (see App C for details of the interpretation of J) . A general 
class of separable potentials consists of linear combinations of 



so that 

V K(3) oc n a - p /{l + K 2 }%+ 2 and V (aj/3) oc sin a - /3 f cos^f . (52) 

The most physically relevant subcase therein are the power-law-mimicking potentials (3 = (remember that I turns out to 
be constant in similarity relational particle mechanics after variation is done) . These correspond to potential contributions 
solely between particles 2, 3. The case (3 = a are potential contributions solely between particle 1 and the centre of mass 
of particles 2, 3 (which is less widely physically meaningful). It also turns out that the duality map sends V( a .f3) to 

V( a ,a-/3)- 

Simple examples therein of quantum-mechanical interest are the constant potential (a = 0) and similarity relational 
particle mechanics' mimicker of harmonic oscillators (a = 2). I choose these for explicit study due to the ubiquity of 
constant and harmonic oscillator potentials in theoretical physics; moreover harmonic oscillator potentials confer nice 
boundedness features at the quantum level. I use the following notation for the constants of proportionality for the single 
harmonic oscillator between particles 2, 3: 2\/^ ^ = h 23 {q 2 — <7 3 } 2 = p\ for h 23 the ordinary position space Hooke's 
coefficient and H\ the relative configuration space's 'Jacobi-Hooke' coefficient. I use the obvious cyclic permutations of 
this for harmonic oscillators between particles 3, 1, and particles 1, 2. Moreover, the linear combination of these last two 
such that rri2hi 3 — m 3 h\ 2 (resultant force of the second and third 'springs' pointing along the line joining the centre of 
mass of particles 2, 3 and the position of particle 1) is 2{v| 2 ^ + V^} = H[ 2 ^ p\ + H^p p\, where 

H[ 2) = {h 13 m 2 2 + h 12 ml}/{m 2 + m 3 } 2 and H ( 2 } = h 12 + h 13 . (53) 

Additionally, one can consider V = V^ 1 - 1 + V^, whereupon one has ffi = h[ 1} + H {2) and H 2 = H {2) ; then define 
Ki = Hi/p,i so that 

2V = H lP 2 + H 2 p 2 = K x i\ + K 2 i\ . (54) 

I refer to this case as the 'special' case of the 3 harmonic oscillator problem, the general case having angle-dependent 
potentials (see the next SSec) due to not having the resultant force of the second and third 'springs' point along l 2 ). 
Three notes of importance below and paper II are as follows. 

1) In (1Z, $) coordinates this leads to the form 

2\/ ^{KiK 2 +K 2 }/{l+K 2 } 3 , (55) 
or, in (6, $) coordinates, by the simplifications cos 2 ^ = {1 + cos0}/2 and sin 2 ^ = {1 — cos0}/2, to the form 

V = A + Bcos9 for A = {K 1 + K 2 } /16 , B = {K 2 - /16 (56) 

2) This situation is the linear combination of a V( 2 .o) an d a V( 2 .2) and as such is self-dual provided that the values of K\ 
and K 2 are also interchanged (or, equivalently, the sign of B is reversed). 

3) if Ki = K 2 (corresponding to the highly symmetric balance m\h 23 = m 2 h 3 i — m 3 /i 12 in position space), then 

V = A , (57) 
to which I refer as the 'very special case'; this is unconditionally self-dual. 

4.3 Scalefree triangleland with harmonic oscillator type potential 

In all other multiple harmonic oscillator cases, there is a ^-dependent cross-term, rendering the potential nonseparable 
in these coordinates. The general triple harmonic oscillator like dipotential 2V = h 23 {q 2 — <Z 3 } 2 + cycles maps to 

2V = {K{R 2 + LTlco^ + K 2 }/{1 + K 2 } 3 for L = {pip 2 y 1/2 2{h 13 m 2 - h 12 m 3 }/{m 2 + m 3 } . (58) 
Or, alternatively, in the spherical presentation, the potential is 

V = A + Bcos9 + Csin6cos$ for C = L/16 . (59) 
This clearly includes the previous SSec's special case by setting C — (and so L = 0, corresponding to 7712/113 = TO3/112). 
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4.4 Analogy with some linear rigid rotor set-ups 

Some useful mathematical analogies for scalefree triangleland with multiple harmonic oscillator like potentials are as 
follows. 

very special harmonic oscillator < — ► linear rigid rotor , (60) 

special harmonic oscillator < — ► linear rigid rotor in a background homogeneous electric field in the axial ('z')-direction , 

(61) 

general harmonic oscillator < — ► linear rigid rotor in a background homogeneous electric field in an arbitrary direction . 

(62) 

In particular, this classical problem has 

T rot or = Irotor{^ 2 + sin 2 0</> 2 }/2 , V rotor = -M£cos9 (63) 

where I ro tor is the single nontrivial value of the moment of inertia of the linear rigid rotor, £ is a constant external electric 
field in the axial 'z' direction and M is the dipole moment component in that direction. Thus the correspondence is 
6 < — ► 6, $ < — ► <f>, 

1 < ► Irotor , (64) 

(energy)/4 - (sum of mass- weighted Jacobi-Hooke coefficients)/16 = E/4 — A = E — A < — > E = (energy) , (65) 

(difference of mass-weighted Jacobi-Hooke coefficients) = B < — ► — A4£ , (66) 

These all being well-studied at the quantum level [691.170]. this identification is of considerable value in solving the relational 
problem in hand, by the string of techniques in Fig 3. 

[Figure 3 In this paper I perform coordinate transformations to obtain standardly solvable problems and then back-track to see what happens 
in terms of the original mechanically-significant variables.] 

4.5 A brief study of the potential 

Working in spherical coordinates, set = <9V/<99 = — BsinO + Ccos0cos$, = dV/d$> — — Csin8sin<i> to find the critical 
points. These are at (0,$) = (arctan(C/S), 0), (— arctan(C/-B), ir) which are antipodal (see Fig 4); in fact the potential 
is axisymmetric about the axis these lie on. The critical points are, respectively, a maximum and a minimum. [The very 
special case B = C = is also critical, for all angles - this case ceases to have a preferred axis.] 

[Figure 4] 

1Z small corresponds to small, for which 

U + E = {E-A-B}- C0cos$ + BQ 2 /2 + O(0 3 ) , or (67) 

2{U + E} = 2E — K 2 - LTZcos^ - {4E + Ki- 3K 2 }K 2 + 0{K 3 ) = Qo - LTZcos^ - {Q 2 1Z 2 } + 0(1Z 3 ) . (68) 

Thus the leading term is a constant, unless Qo = 2E — K 2 ( cx E — A — B ) = 0, in which case it is linear in or 1Z 
and with a cos$ factor, unless also L ( cx C ) = (which is also the condition for the 'special' case), in which case it 
is quadratic in O or 1Z, unless B — (given previous conditions, this is equivalent to Q 2 = 4E + K\ — 3K 2 — 0), which 
means that one is in the K 2 = 2E subcase of the 'very special' case, for which U + E has no terms at all. 
1Z large corresponds to the supplementary angle 3 = tt — Q being small, so 

U + E = {E - A - B} + CScos$ + BE 2 /2 + 0(E 3 ) , or (69) 

2{U + E} = {2E-K 1 }TZ- 4 -LTZ- 5 cos^-{4E+K 2 -3K 1 }TZ- 6 +0{TZ- 7 ) = Q 4 7^ 4 -£7^ 5 cos$-Q 6 7?.- 6 + C>(7^ 7 ) . (70) 

Thus the leading term goes as a constant in S or as TZ~ 4 , unless Q4 = 2E — K\ ( cxE — A + B ) = 0, in which case it goes 
linearly in 5 or as 1Z~ 5 in each case also with a cos$ factor, unless also L ( cx C ) = ('special' case), in which case it 
goes quadratically in 2 or as 1Z~ 6 , unless B — (given previous conditions, this is equivalent to Qq = 4E + K 2 — iKi = 0), 
which means that one is in the K\ = 2E subcase of the 'very special' case, for which E + U has no terms at all. Note that 
B = implies K\ = K 2 , so this very special subcase indeed coincides with the previous paragraph's. Finally note that 
the large and small asymptotics are dual to each other (the difference of 4 powers is accounted for by how the kinetic 
energy scales under the duality map) , so that one need only analyse the parameter space for one of the two regimes and 
then obtain everything about the other regime by simple transcription. 
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4.6 Classical equations of motion 

The Jacobi-type action for this problem is, in spherical coordinates and using the barred banal conformal representation, 

f 1 9 2 + sin 2 6$ 2 f E 1 
S = / dXd —- I - - A - Bcos9 - Csin8cos$ V . (71) 

Then the Euler-Lagrange equations are 

9^ - sin9cos9{$*} 2 = 5sin9 - Ccos9cos$ , (72) 

{sin 2 8$*}* = Csin6sin$ . (73) 
On account of the action being independent of t, one of these can be replaced by the energy integral 

{9^} 2 + sin 2 e{^} 2 +A + BcQse + CsinQcos<J) = E (74) 

A dynamical systems and phase space analysis of these equations and their interpretation in terms of Ii, I2 and <!> will be 
presented elsewhere [71]. Also note that the notions of V e g- and V or b are inapplicable in the nonseparable case. 

4.7 'Special case' 

For C = 0, d reads 

If- 7 2 ~) - 

'{e*} 2 + ^2— }+A + Bcose = E (75) 



2 [ sin~9 
which can be integrated [using (j43| to eliminate i cm in favour of $]: 

$ - $0 = J Jde/sinG^E -A- 5cos9}sin 2 9 - J 2 . (76) 

For the very special case B = C = 0, this reduces to the well-known integral for the computation of geodesies for the 
sphere - c.f. [17] , which exact solution I now cast below in terms of 7Z, $ and then in terms of the 'original' t l7 l 2 variables. 
I've also solved the B 7^ case exactly by e.g. tan^ = 7Z = \Jx — 1, obtaining then via Maple [72] a composition of 
polynomials, roots and elliptic functions in x that I consider to be too complicated to present here. 

[T7] sketches V e ff and U or b for the single harmonic oscillator and their large and small asymptotic behaviours. J = 
is rather simpler: 1-d motion i.e. the orbits are straight lines. Assume J ^ from now on in looking for orbits with 
less trivial shapes. The universal large 1Z asymptotic solution's V e ff exhibits a finite potential barrier and the orbits 
are bounded from above. The small 1Z asymptotic solution for a fixed E > harmonic oscillator problem is the usual 
radial/isotropic harmonic oscillator problem. Here, V c ff is an infinite well formed by the harmonic oscillator's parabolic 
potential on the outside and the infinite centrifugal barrier on the inside, and the orbits are bounded from below. For 
a fixed E > harmonic oscillator scalefree triangleland problem, in contrast with the usual radial harmonic oscillator 
problem, the potential tends to rather than +00 as 1Z — > 00, and the orbits are bounded from both above and below. 

4.8 Exact solution for the 'very special case' 

This effectively constant- V case is equivalent to the geodesic problem on the sphere. Thus it is solved in 9, $ variables 
by great circles (first equality) 

J Fcos($ - $ ) = 2cot9 = {1 - 1Z 2 }/1Z = {i\ - il}/LiL 2 (77) 

(following through with variable transformations so as to cast it in terms of straightforward relational variables, and using 

JF = 2\J 2{E — A] - J 2 = V2E -K 2 - 4J 2 = ^Q Q - AJ 2 ^j . The 1Z form above can be rearranged to the equation of 

a generally off-centre circle [Fig 5a)]. 

While, in terms of the original quantities of the problem, 

F 2 {Li ■ L2V + I kill 4 + lk 2 H 4 + 2F{|y| 2 - |U 2 || 2 Hi ^ 2 cos$o = {2 + F 2 sin 2 <I> }|y | 2 |y | 2 . (78) 

Note 1) this solution is homogeneous in ||tx||, IU2II an( i \JL\ ' L2 ( m f ac t it is a fourth-order homogeneous polynomial 
in these quantities). These quantities are all invariant under rotation, so the zero total angular momentum constraint 
is manifestly satisfied, while the zero dilational momentum constraint is manifestly satisfied due to the homogeneity. 
['Homogeneous equation' is another way of expressing 'equation in terms of ratios alone'.] 
2) It has simpler subcases: sin"I>o = gives the merely second-order 

^i-i 2 + IUill 2 = IU 2 l| 2 . (79) 
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which I sketch in Fig 5, while cos$o = gives 

F 2 W-L 2 r + I W I 4 + IU2II 4 = {2 + ^IIUilHy I 2 • (80) 
3) The solution is invariant under tj < — > t 2 provided that also $0 is shifted by a multiple of n. 

[Figure 5 a) On the sphere the solutions are great circles, while their projections onto the (TZ, <!>) plane are circles. I also provide a sketch of 
how the projection works; the equator maps to the circle TZ = 1 and the 'Greenwich meridian' maps to the vertical line. 

b) The barycentric partial moments of inertia as functions of $ are then [by H77II . that their sum is a constant I, and that this sketch is 
readily rotatable from <&o = to <&o arbitrary and using [72] ] Ii , I2 = i |l± Fcos«I>/v / 4 + F 2 cos 2 $|. This is additionally, as expected, readily 
rescaleable for any I ^ 0, so I take 1 = 1. For diagrams of this kind in this paper, I use the convention of black for Ii and dashed for I2. Then 

as F > 0, one gets twice the circle of radius 0.5 (corresponding to the TZ = 1 solution); the first picture, for F = 0.1, is well on its way toward 

that behaviour. As F grows, there motion has two 'almost- halves' in which one of Ii, I2 completely dominates over the other, separated by 
two thin wedges in which Ii, I2 cross over. This amounts to I1T2 undergoing periodic oscillations. Note that they always cross over at 1/2, 
corresponding to all the other circles in Fig 5a) cutting the TZ = 1 circle. Also note the complete symmetry between the status of Ii and I2 

in the opposite direction (<E> > <E> + 7r). As F ► 00, the two curves touch the origin and, oppositely, the unit circle, corresponding to TZ 

running from to 00, i.e. the vertical line of Fig 5a). 

c) Re-intcrprctation in terms of particles. I consider only mi = = for simplicity. The equator (TZ = 1 circle) F = case corresponds to 
going through the configurations sketched, which can be summarized by the picture underneath particles 2, 3 fixed and the circle representing 
the range of positions that particle 1 takes. I also provide the analogous summarizing sketch for the slightly inclined great circle of Fig 5 
a) (this is for $0 = 0, in general different values of t&o correspond to different particle behaviours). The Greenwich meridian corresponds to 
motions between the double collision of particles 2, 3 and the 2, 1, 3 collinearity with particle 1 at the centre of mass of particles 2, 3. 

d) Some useful terminology for the particle configurations. Some definitions to aid this are as follows: near- isosceles and near-collinear (for 
a < some qq, constant), ' Jacobi-flat' (Ii << I2), ' Jacobi-regular' (Ii = I2, which, when isosceles, is equilateral) and 'Jacobi-tall' (Ii >> I2).] 



4.9 Small relative scale asymptotic behaviour 

The first approximation for TZ is the generic asymptotics (Qo = 2E — K 2 ocE — A + B ^ 0, soU + E goes as Qq/2). [Note 
however that not all dynamical orbits enter such a regime - sometimes the quadrature's integral goes complex before the 
small TZ regime is attained (TZ large 'classically forbidden') when \J(TZ)TZ 2 — J 2 < 0.] Then integrating the TZ quadrature 
([50l) gives the orbits 

± V2^sec($ - l>) = TZ = tanf = l x /l 2 (81) 



(following through with variable transformations so as to cast it in terms of straightforward relational variables, and using 
the J-absorbing constant q = Qo/2J 2 = E - K 2 /2/J 2 = 4{E -A- B}/J 2 ). Note that the TZ form of the orbits are 
parallel straight lines (vertical for $0 = 0, 7r) [Fig 6a)]. Moreover, these are known not to be very good approximands in 
that they totally neglect the non-constant part of the potential and thus are precisely rectilinear motions. Thus at least 
here (and one may suspect likewise in the next section), it is often necessary to use the second approximation in studies. 
The other equalities in (|8ip convert that result into what form the asymptotic orbits take in straightforward relational 
variables; I sketch these in Figs 6a) and 7a). 

While, in terms of the problem's original quantities, 

||t 2 || 4 + 2g { a -^} 2 = 2q {|y| 2 ||^|| 2 sin 2 $ + 2y2^||i 2 || 2 {ii -t2} 2 }cos$o , (82) 

which is again a fourth-order homogeneous polynomial in \\l 2 \\ an d ^Jk\ ' L 2 ■ It admits the simpler cases 1) sin$o = 

(which is merely second-order): 

= V^QoLi ■ ki , (83) 



||t 2 || = V^?o£i • £2 > 
and 2) cos$ = 0: 

IU2H 4 = 2go{|U 1 || 2 |U 2 || 2 -{£ 1 ^ 2 } 2 }. (84) 

The second approximation for 1Z small in which L = but the Q 2 TZ 2 /2 term is also kept turns out to also often be 
necessary. If Q = but L^0 one has non-generic asymptotics not directly covered in this SSec, but one can use the 
technique of Sec 5 to reduce this case to one of those covered in the present SSec in a new set of variables. If Qo, L = 0, 
one resides within the special case, and the next leading term is Q 2 1Z 2 /2, which case is included in my working below; 
indeed now the "second" approximation is always necessary. If Q 2 = 4E + K\ — 3K 2 = 16{E — A— 2B} — also, one is 
within the very special case, and so one does not need any asymptotic calculations as one has the exact solution of SSec 
4.8. 

For the second approximation integrating the 1Z quadrature gives the orbits 



± 1/y q + ^q 2 - g 2 cos(2{$ -$})=■£ = tanf = n/i 2 (85) 

(following through with variable transformations so as to cast it in terms of straightforward relational variables, and using 
the ^/-absorbing constant q 2 — Q 2 / J 2 ). This is straightforwardly rearrangeable into quite a standard form (e.g [73l 174) ) 

TZ 2 = 1 , (86) 

lo + V3o - 92Cos(2{$ - $0}) 
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the case- by-case analysis of which is provided in Fig 6b) . 

While, in terms of the problem's original quantities, and using g for \/<7 2 — q 2 , 

\\L 2 \\ S + Uo +3 2 cos 2 (2$ )}|M 4 |y | 4 - 2 <70 |U 2 || 6 |y | 2 + 4 3 2 { a - ^} 4 + 

2 5 |U 2 || 2 {9o||^i|| 2 - || i 2|| 2 }{2U 1 - A2 > 2 - |U 1 || 2 |U 2 || 2 }co S (2ci> ) ^ ^g 2 ^-^} 2 !!^!! 2 !!^!! 2 , (87) 

which is an eighth-order homogeneous polynomial in H^H, ||6 2 || and • i 2 . It admits the simpler cases 1) sin(2$o) = 
(which is merely fourth-order): 

IU 2 || 4 = te-<?}lkiH 2 |U 2 || 2 + 2 5 U 1 -, 2 } 2 (88) 

and 2) cos(2$ ) = 0: 

IU 2 || 8 + go 2 lkill 4 lk 2 ll 4 -2 g olk 2 || 6 lk 1 || 2 = V^.^iii^ini^H 2 -^.^} 2 }. (89) 

Taking qo = further simplifies the simple subcases. While, taking 5 = collapses the general solution to 

||t 2 || = (2, 1, 3 collinearity with 1 at the centre of mass of 2, 3) or||t 2 || = iy^golli-JI (rectilinear motion: $ fixed) . 

(90) 

[Figure 6. a) The vertical straight lines of the first small approximation. 

b) qo > 0, <?2 > (which is the usual mathematics of the isotropic harmonic oscillator), gives ellipses centred about the origin 175) . up to 
maximum values on the 52 = <?g parabola on which they become circles centred on the origin; beyond this the solution ceases to exist. While, 
go > 0, q2 = recovers pairs of straight lines belonging to a). 92 < is the upside-down isotropic harmonic oscillator and gives hyperbolae, 
which are obtuse for qo > rectangular on 50 = an d acute for go < 0- Finally, qo = ?2 = is the circle at infinity, while go < 0, g2 > and 
go = 0, 52 > are classically forbidden, as can be seen from the quadrature having to be real. 

c) The partition of parameter space into these cases; it is obvious from (50) that the shaded region is classically disallowed.] 

[Figure 7. Using a form valid for both first and second small approximations, the partial moments of inertia are, for $0 = 0, 
Ii = I/{1 + go + gcos(2$)} and I 2 = I{g + gcos(2<t>)}/{l + g + gcos(2$)}. 

a) Then the three types of behaviour of the first large approximation are: 'I2 is surrounded by Ii', 'Ii and I2 touching', and 'Ii and I2 cross-over' 
(of which I provide two representatives). The touching case is the limiting case between the other two cases. The only instance in which the 
small approximation is self-consistent is in the last picture for <!> within wedges around five angular tick-marks wide either side of and of it. 

b) Next, here are sketches of the nine further behaviours exhibited by the second approximation. The first three are limiting behaviours on the 
g2 = g^ parabola, of which the second is itself the limiting behaviour between the other two a single point (corresponding again to the 1Z = 1 
solution). The next three are the unfolded version of 'cross-over', the reverse case of 'touching' (again a limiting behaviour) and the reverse 
case of 'is surrounded'. The last four are all cases in which there occurs 2, 1, 3 collinearity with particle 1 at the centre of mass of particles 2, 
3: 'is surrounded', a 'touching' limiting case, and two instances of 'cross-over'. The outwards-lying arrows indicate for which 3> in each case 
the small approximation is self-consistent. 

c) Here I provide a sketch (not to scale) of which regions of the parameter space these various cases reside in. The dashed curve touches the 
parabola, thus trisecting the classically-allowed region into: a region above and to the left of the dashed curve where the small asymptotics is 
everywhere-valid, a region to the right of the dashed curve where it is valid in some wedges, and a region below and to the left of the dashed 
curve where it is nowhere- valid. 

d) A few examples of reading off what is happening in the particle position picture. The third and sixth cases in 7b) correspond to particle 
1 describing a closed curve round particles 2, 3. The fifth case describes particle 1 coming in from infinity (= particles 2, 3 colliding) and 
the approximation breaking down for some $ c . The tenth case describes particle 1 coming in from infinity and reaching the centre of mass of 
particles 2, 3 at some angle <& c .] 



4.10 Large relative scale asymptotic behaviour 

For analogous notions of first and second large approximations, now the quadrature in U = 1/1Z takes the same form as 
the 7?.-quadrature in the above workings with 

Qo — > Qa , Q2 — ► Q 6 (91) 

(94, 96 and / below are then the obvious analogues of go, 92 and g). Hence the solutions are dual to those of SSec 4.9's. 
Thus all of SSec 4.9's results apply again under the duality substitutions (and the subsequently-induced language changes 
'small' — > 'large', and "2, 1, 3 collinearity with particle 1 at the centre of mass of particles 2, 3" — ► "collision between 
particles 2, 3, also interpretable as particle 1 escaping to infinity"). In particular, the T,I 2 plots of Fig 7b) and which 
parts of the parameter space the various cases hold in and where the now small asymptotics is valid can just be read off 
the existing figures under these substitutions, and then the region in which small asymptotics applied before [Fig 7c)] is 
now that in which large asymptotics now applies. 

New 1Z, $ plots are, however, required (Fig 8) Some particular comments are: 1) the first approximation is then 

± v /2<fccos($ - $ ) = Tl = tanf = l 1 /l 2 ■ (92) 

In the (71, $) plane and for <!>o = 0, this takes the form of a family of circles of radius ^Jq^/2 and centre (1/54/2,0), so 
that they are all tangent to the vertical axis through the origin ([T7] and Fig 6a)). 2) The second approximation gives 



±d qi + Jql- g 6 cos(2{$ - *}) = Tl = tanf = Li /l 2 . (93) 
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In the 1Z, 4> plane and for <£>o = 0, this take the forms in Fig 8b) in the parameter regions delineated by Fig 6c). 

[Figure 8 a) The first large approximation's family of tangent circles in the (Ti, $) plane. 

b) Additional behaviours of the second large approximation: bulging, rectangular and thin tear drops, ellipse-like and peanut-like curves and 
circles centred on the origin. The partition of qi, qe parameter space here is the same as that of goi <?2 parameter space in Fig 6c); indeed Fig 
6c) is an important clarification of the behaviour exhibited by the large approximation.] 

Note: for general V( Q ) with constant of proportionality A Q , one gets exactly the same large-asymptotics analysis as 
here, with qo = 2E — A a , <?2 = 4E — {4 + a}A Q . Thus generally this duality to the isotropic harmonic oscillator of the 
universal large-scale asymptotics of scalefree triangleland is a useful and important result for this Machian mechanics, for 
the classical and quantum mechanics of isotropic harmonic oscillators is rather well-studied and thus a ready source of 
classical and quantum methods, results, and insights. 



4.11 Investigation of the intermediate-7^ region 

Note that only two of Qo, Q2, Qi and Qq are independent. Thus prescribing a particular small asymptotics entails 
prescribing a particular large asymptotics too. There are however no compatibility restrictions: each small asymptotic 
behaviour is capable of connecting to each large asymptotic behaviour. 

Also note from the twin circles, tear drops and peanuts that more than one region in which large asymptotics applies is 
possible. Then numerical integration using Maple's [72] rkf45 solver reveals that precession can occur in the intermediate- 
1Z region, so that orbits can have arbitrarily many large (and small) asymptotic regions. However, the spiraling reported 
in [T7] does not match this paper's asymptotic solutions in all relevant places. I account for this as follows. 1) I observe 
that when the bona fide energy constant E ^ 0, solutions do include the circle centred on the origin and spirals to/from 
that [73] which resemble the previously-reported spiraling. 2) However the scalefree triangleland problem in hand has E 
= and so possesses no such solutions. 3) The previously-reported spiraling is due to |17j's code not implementing the 
E = restriction. Thus the observation of spiraling with that code is unsurprising, and wrong. The code I used for the 
present paper specifically builds its 7?.-derivative 'initial datum' by solving the E = energy equation for the input 1Z 
initial datum (and checks that the evolution does not take one away from E = 0), thus amending the error. 

5 Normal coordinates for scalefree triangleland with multi-harmonic os- 
cillator type potentials 

5.1 A rotation sending the general case to the special case in new coordinates 

One can avoid having a C-term by performing a rotation or normal coordinates construction, which preserves the form of 
the kinetic term while sending the potential to a C-free form in the new coordinates. One can get to these by inserting 
such a rotation between Step 1 and Step 2 in Sec 3.2, or by rotating coordinates at the level of the equations of motion. 
I denote the new coordinates with N-subscripts (N for normal). Due to the general case in N-subscripted coordinates 
taking the same form as the special case in the original coordinates, we can uplift Sec 4 to an exact solution of general 
multiple harmonic oscillator case by inserting N-subcripts and, at the end of the calculation, rotating back from normal 
Jacobi coordinates to a more complicated form in terms of the original Jacobi coordinates. 

What is the requisite rotation angle? From the matrix equation i?(through angle a about y axis)n = n N , 

/ cosa -sina \ / C \ / \ 
I 010 I 1 J / \ 

I sina cosa J I B J = \ B N J , (94) 

where n, n N are unit vectors considered to be in unit-radius spherical coordinate form about the original and normal- 
coordinate axes, the requisite rotation is by 

a = arctan(C*/ J B) : (95) 
this sends BcosO + CsinOcos^ to -BncosOn- Then 



sina = C/y/B 2 + C 2 , cosa = B/y / B 2 + C 2 , B N = V ' B 2 + C 2 , (96) 



cos9 N = {BcosO + Csinecos^j/v/S 2 + C 2 and $ N = arctan U/B 2 + C 2 sin0sin$/{Bsinecos$ - Ccos6}J . (97) 
The potential is now 

{Kfp 2 lN + K$f$J/8l = A N + BncosOn . (98) 
It is also useful to note for later use the following coefficient interconversions: 



A N = A , ,K? = 8{A - ^B 2 + C 2 } , X 2 N = 8{A + V ' B 2 + C 2 } , (99) 
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Kf = Ki + K 2 - y/ild-Ktf + L 2 } , = K 1 +K 2 + \K\ — K 2 } 2 + L 2 } . (100) 

From the spherical perspective, the normal coordinates solution has the same form as the special solution in the original 
coordinates, but now one is to project onto the general tangent plane rather than the tangent plane at the North Pole, 
interpreting now that general stereographic coordinate as the ratio of the square roots of the barycentric partial moments 
of inertia. This permits graphical sketches of the qualitative behaviour rather than fairly lengthy analytical expressions 
constructed by passing from (On, $n) coordinates to (©, $) coordinates and then on to mechanically significant variables. 

5.2 Examples: preamble. 

Note that the very special case's solution is invariant under this rotation: that has no preferred axis. So one needs to look 
slightly further to obtain nontrivial examples. First I give the simplest nontrivial example of the first small asymptotics 
from both an analytical and a graphical presentation. For this paper to be of manageable length I then only provide a 
graphical perspective for the first large asymptotics, and brief comments on further cases. 

[Figure 9 N-large and N-small domains of applicability map to the original (TZ, $) plane as indicated.] 

5.3 First small asymptotics solution for general case 

From (|5Tj) with N-subscripts, tan^p = sin9/{ 1 + cos6n} (|97 ll99|l and elementary cancellations, for ^ = the analytic 
solution for this takes the form 

B N + Scos© + Csin6cos$ = y / 2^°"{Bsinecos$ - Ccos©} , (101) 

where bears the same relation to K 2 as qo bears to K 2 . Then in terms of (1Z, $), 

B N -B- y^C* j TZ 2 + 2 jc - y^sj ftcos$ + B N + B + y^c = . (102) 

Or, in terms of straightforward relational variables T, l 2l $, 

B N - B - y^cj T + js N + B + y^cj I 2 + 2 jc - \/2^s| ^/hhcos® = . (103) 

Finally in terms of the original variables for the problem, 

B N -B- y^cj \\ Ll \\ 2 + js N + B + yS^cJ |y | 2 + 2 jc- y^sL -i 2 = . (104) 

Like for the C = case (|83|) , this is a second-order homogeneous polynomial in 1 1 l x \ | , 1 1 1 2 1 1 and ■ i 2 . 

[Figure 10 a) Sketch of how the first small approximation's parallel straight lines in the (R., ©) plane project onto the (©, $) sphere. 

b) (7?., "3?) sketches for C ^ 0, either from plotting the analytical function or from projecting the sketch on the sphere onto the tangent plane 
of the appropriately-rotated North Pole. The family of parallel straight lines for C = is now a family of circles tangent to a single point. 
Within the (unshaded) region where this solution is valid, one therefore obtains circular arcs. 

c) Subsequent sketches of Ii and I2 as functions of 3> showing some of the distortions which occur when a C-term is switched on in the 
'cross-over' case of Fig 7a). All subfigures have A = 1 = B and E = 9. The first picture is for C = 0.1. The second picture is for C = 0.5, 
by which stage I2 encloses the origin. The third picture is for C ~ 2.9, for which Ii and I2 touch, after which I2 'is surrounded' by Ii (fourth 
picture).] 



5.4 First large asymptotics solution for general case 

[Figure 11 a) Sketch of how the first large approximation's family of tangent circles in the (Ti, 0) plane project onto the (©, <t>) sphere, 
b) (Ii, <J>) sketches for C ^ either from plotting the analytical function or from projecting the sketch on the sphere onto the tangent plane 
of the appropriately-rotated North Pole. One still has arcs that touch at a point, but this point is no longer at the origin but is rather shifted 
away from it along the $ = axis. This does not cause any noteable changes to Ii and I2 as functions of <E>, so I do not provide sketches of 
these.] 



5.5 Further examples 

One can go on from here to provide in rotated coordinates general solutions to the second small and large asymptotics, 
and of the general numerical behaviour in the region in between. Here are some brief comments. 

1) Similarly to the C = approximation having turn-around ellipses rather than having to follow straight lines, the second 
small approximation allows for turn-around behaviour rather than having to complete the first small approximation's 
circles: approximate circular arc, turn-around, approximate circular arc in opposite direction to the original. 

2) One can now have now asymmetric bulges where one had symmetric ones before (e.g. in the first large approximation 
or in the peanut case of the second large approximation) ; indeed one bulge can non-generically become infinitely big (like 
the straight line in the first large approximation) and even form another contribution 'beyond infinity' which shows up 
on the other side of the opposite bulge. 
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6 Scaled triangleland at the classical level 

6.1 Scale— shape coordinates (1,7^, $) or (I, 0, $) 

The general scaled triangleland with multi-harmonic oscillator like potential's Jacobi-type action is 



s| u-,,,,-,.„ , , ,, VtI e + u, s W w } 



The Euler-Lagrange equations are 



*2 , „:„2q**2i 1 f 5V , E + U \ _n ^0*1* T 3**3..- ~ ~ 1 ^ 



r* - i{e*" + sin"e$*"} + — |— + —j— | = o , {re*}* - r$* 2 smBcosB + — — = o , (106) 

{I 2 sinWr + l^ = 0, (107) 

and there is an accompanying energy integral 

{I* 2 + I 2 {B* 2 + sm 2 B<T 2 }}/2 + V/4I = E/4I . (108) 

[Above, 



In the case in which V is independent of $, (|107[) simplifies to the first integral 

I 2 sin 2 B$* = J . (110) 

See App C for physical interpretations of the relative angular momentum quantity J. 
Then the other Euler-Lagrange equations and the energy integral take the forms 

OS j2 1 P V E - V \_n fl 8r*" J 2 C0S6 . 19V 









^ + 


I* 2 




~2~ 


+ 2 



i**-ib* 2 -^ — =- + — <^ — + — — } = o , {i^en* --5 — 5- + tf^ = o, (in) 

I 3 sin 2 6 41 I dl I J 1 1 I 2 sin 3 9 41 99 ' v ' 

J 2 V _ JE 

' 2I 2 sin 2 6 + 41 ~ 41 ■ (112) 

6.2 Scaled triangleland with general multi-harmonic oscillator potential 

The usual general multi-harmonic oscillator potential maps to 

{l/4I}I{{iTi + K 2 }/4 + {K 2 - K^/^cosQ + {L/4}sinBcos$} = A + ScosB + CsinBcosB . (113) 

6.3 Sketch of the potential 

The sketch of potential is as for the corresponding similarity problem. The sketch of V — E is slightly different from that 
of V — E for the similarity problem in that one has a further equation in looking for critical points, 

<9{V-E}/<9I= E/4I 2 , (114) 

so one also needs E = in order to have these. After that the analysis is as before (including the picking out of a preferred 
axis except in the B = C = case), except that the Hessian has an extra row and column of zeros. 

6.4 Classical equations of motion 

The Jacobi-type action for this problem is then 



S = 2 / dA^ 2+l2 ^ 2 2 +SiD2e<t)2 > { I - A - BcosB - CsinBcos*} . (115) 
The Euler-Lagrange equations are 

r*-i{e* 2 + sin 2 e$* 2 } + E/i = o , (ne) 

{I 2 6*}* - I 2 $* 2 smBcos6 + C*cos6cos$ - Ssin9 = , (117) 

{I 2 sin 2 e$*}* - CsinOsin^ = . (118) 

The energy integral is 

I* 2 + I 2 {6* 2 + sin 2 9$ i2 } E , . 
i +A + BcosG + Csinecos$ = y . (119) 

For small 1Z or B and for large 1Z or small 5 = ir — B, Sec 4's results for approximations to the potential hold again. 
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6.5 Special Case 

In the special case of C = 0, p07p applies and the remaining Euler-Lagrange equations and energy integral become 

I** _ ie*2 _ j 2 /i 3 s in 2 e + E/I = , {I 2 6*}* - J 2 cose/I 2 sin 3 e + C*cos9cos$ - BsinG = , (120) 
I* 2 /2 + I 2 9* 2 /2 + J 2 /2I 2 sin 2 6 + A + BcosO + Csin6cos$ = E/I . (121) 

6.6 Analogy between very special case and Kepler— Coulomb problem 

The very special Euclidean relational particle mechanics harmonic oscillator banal-conformally maps to the Kepler problem 
with 

(radius) = r < — ► I (total moment of inertia) , (122) 
(test mass) = m < — > 1 , (123) 
(angular momentum) = L < — ► J (relative angular momentum - see App C) , (124) 
(total energy) = E < — > — A = - (sum of mass- weighted Jacobi-Hooke coefficients) /16 (125) 

and 

(Newton's gravitational constant) (massive mass) (test mass) = GMm < — > E (total energy)/4 (126) 
[or to the 1-electron atom Coulomb problem with this last analogy replaced by 

(nuclear charge) (test charge of electron)/47r(permettivity of free space) = (Ze)e/47reo < — ► E (total energy)/4 ] . 

(127) 

Also note that the positivity of the Hooke's coefficients translates to the requirement that the gravitational or atomic 
energy be negative, i.e. to bound states. While, the positivity of E required for classical consistency corresponds to 
attractive problems like the Kepler problem or the atomic problem being picked out, as opposed to repulsive Coulomb 
problems. 

Also, the special case corresponds to the same 'background electric field' that the rotor was subjected to in Sec 4, 
which, moreover, is proportional to cosG which is analogous to cos#, which is in the axial ('£') direction but is not 
the well-known mathematics of the axial ('z') direction Stark effect for the atom, which involves, rather, rcosO. But, 
nevertheless, the situation in hand is both closely related to the rotor situation in Sec 4 and to the mathematics of the 
atom in parabolic coordinates (see e.g. [751(70]). The general case is then the same situation but with the 'electric field' 
pointing in an arbitrary direction. The idea then is to use the obvious analogue of the scheme in Fig 3 to solve Euclidean 
relational particle mechanics problems in straightforward relational, relative and absolute terms. 

6.7 Exact solution for the very special case 

Now 

I i2 /2 + J 2 /2I 2 sin 2 9 + A = E/I (128) 

[usually one would set 0o = 7r/2 without loss of generality, however the present physical interpretation has the value of 
0o be meaningful, as ©o = 2arctan(ii/t2)]. Thus the solutions are conic sections 

I = //{l + ecos($-$ )} (129) 
where the semi-latus rectum and the eccentricity are given by 

I = J 2 /Esin 2 e , e = \Jl- 2AJ 2 /E 2 sin 2 e . (130) 
So, in terms of straightforward relational variables, 

Ii+I 2 = //{l + ecos($-$ )} , (131) 
and in terms of the original variables of the problem, 

{Ikill 2 + IN PHIkilliy I + • i 2 cos$ + V\\l 1 \\ 2 \\l 2 \\ 2 -{l 1 -L 2 } 2 ^o}} = I ■ (132) 

Note that this is non-homogeneous (it rearranges to a non- homogeneous eighth order polynomial) ; this is OK as solutions 
of Euclidean relational particle mechanics do not have to be scale- invariant. 

In moment of inertia-relative angle space, for 2A = {Esin0o/J} 2 one has circles, for < 2A < {EsinOo/J} 2 one has 
ellipses, for A = one has parabolae (corresponding to the case with no springs). The hyperbolic solutions (A < 0) are 
not physically relevant here because this could only be attained with negative Hooke's coefficient springs. The circle's 
radius is I = I = E/2A while for the ellipses I is bounded to lie between J/\/2Asin8o and E/2A. The smallest I attained 
in the parabolic case is J 2 /2Esin 2 0o- The period of motion for the circular and elliptic cases is 7rE/\/2A 3 . 
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As regards the individual subsystems, combining the fixed plane equation and the I($) relation, 

Ii = /sin 2 f /{1 + e cos($- $ )} » I2 = ^os 2 f /{l + e cos($ - $ )} (133) 

so each of these behave individually similarly to the total I. In the 80 = 7r/2 plane, they are both equal (and so equal 
to 1/2). Circle and ellipse cases have Ii and I2 as closed bounded curves which sit inside the curve that I traces. The 
parabolic case has Ii, I2 curves to the 'inside' of the parabola that I traces. 

6.8 Special case solved 

Use Ii, I2, $ coordinates given by 

Ii = {I- Z}/2 = I{1 -cos6}/2 and I 2 = {I + Z}/2 = I{1 + cos6}/2 , (134) 

which invert to 

I = Ii + h , 9 = arccos ({I 2 - Ii}/{Ii + I 2 }) (135) 

and are mathematically parabolic coordinates scaled by 1/2, which moreover in the present relational context have the 
physical interpretation of partial moments of inertia of the two subsystems. Then 



S = 2 J dA^T{E+ U} = 2 J dA A 



1 f i\ f 2 4i 1 i 2 $ 2 1 r E K x \ + K 2 l 2 



h Ii + 1 

for T, U, E as before. Then the <E>-Euler-Lagrange equation is 

4IiI 2 $" 



(136) 



I1+I2 



= J , (137) 



and the energy integral is, subsequently, 

If If J 2 f 1 11 K x \ + K 2 l 2 E 



2I X 2I 2 8 \li I 2 J 8 4 ' 

which separates into 



(138) 



4I* 2 + J + Kdf = 2EjIj (139) 

for Ei + E 2 = E. This is solved by 

t - t = {2/^} arccos Ul^ - E. t }/^E 2 - K. t J 2 ^) (140) 

(in agreement with |17j . once differences in convention are taken into account). Thus, synchronizing, one part of the 
equation for the orbits is 



V^ 2 arccos I {1^ - Exj/^/E 2 - K lt P j = V-Kiarccos 1^{I 2 K 2 - E 2 }/^E 2 2 - K 2 ,P j . (141) 

[One can see how the arccosines cancel in the very special case... Then Ei = E 2 = E/2 gives L\ = l 2 i.e. 6 = 

2arctan(ti/t 2 ) = 2, arctan(l) = 7r/2, so arc confined to the plane perpendicular to the chosen Z-axis.] Then the <&- 
Euler-Lagrange equation implies 



* - *o = J fit (1 + i| = f VKl [ ^ = V arctan M - ^ ~ 

J \li h) 2f^ V l J Ficosri + Ei ^ I V {Ei + Fi}{F< + A,} 



(for Ti = 2\fKi{t — <o}, Fj = J Ej — KiJ 2 and = Kili — Ei, ), which simplifies to 



* - $ = Yl arctan U {{\/E- -^J 2 - ^ J h + J2 }/ { { \/ E ?-^ j2 + Et j I 2 - J 2 } j (142) 



in the straightforward relational variables. While, in the original variables of the problem, 



X 2 arccos {|M | 2 ^i - E^/JEf - KiJ 2 = ^arccos {\\l 2 \\ 2 K 2 - E 2 }/a/E 2 - K 2 J 2 , (143) 



2 



arccos 



i=l 



^-rr U$0+E arCtan h/ V E * ~ ^ ^ ~ Ei |2 + ^ / 1 V E ? - ^ P + E * H^ll 2 - ^ 



(144) 
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6.9 The single harmonic oscillator case requires a separate working 

For Ki = K 2 = 0, the trajectories are given by, after some manipulation, 



VE 2 /EiI 2 = Ii = scc(E 1 {$ - $o}/E)/V2ET , 

which is obviously the expected straight line in the absense of forces. 

For K 2 = 0, K\ 7^ 0, the trajectories are given by, in straightforward relational variables, 



{IxKx - Eij/^E? - KlP = cos ^2K 1 /E 2V / 2E 2 l 2 - 



(145) 



(146) 



and 



$ — $o = arctan 




E 2 - A'i J 2 — Ei > Ii + J 2 



E 2 -K 1 J 2 + E 1 \l 1 - J 2 



arctan 



'2E2I2 

J 2 



While, in terms of the original variables of the problem, 



{Kiiy | 2 - EJ/a/E 2 - KiJ 2 = cos [ v /2if 1 /E 2V / 2E 2 ||t 2 || 2 - J 2 

$o + arctan 



IMIIMI 



-2||i2l 
J 2 



arctan 




E 2 — ifi J 2 — Ei > ||i;i|| 2 + J 2 



E 2 1 -K 1 J 2 + E 1 \\\ kl 



J 2 



(149) 



6.10 A brief interpretation of the previous two subsections' examples 

In SSec 6.9's example, l 2 (or I 2 ) makes a good time-standard as the absolute space intuition of it 'moving in a straight 
line' survives well enough to confer monotonicity. It is convenient then to rewrite (|1461fl47|) as a curve in parametric form 
with I 2 playing the role of parameter, leading to the plots in Fig 12. 

[Figure 12 a) 3-d plot showing oscillatory behaviour including some changes in the size of the relative angle that occur at regular intervals 
but can involve 'sporadic' changes in how much the relative angle changes in each interval. (The particular plot given is for the 1 harmonic 
oscillator case, with K\ = Ei = E2 = 1 and J = 0.1, with $ plotted vertically, I2 out of the page and Ii into the page.) This is because 
although 1 is 'moving in a straight line' in absolute space, the position with respect to which its separation is being measured from in relational 
space (the centre of mass of particles 2, 3) is then also moving around due to the oscillations of the 'spring' between these particles, 
b) Polar plot of Ii and I2 as functions of <t> for the first oscillation. Further oscillations correspond to similar angular variations at larger radius. 
N.B. that Ii and I2 are independent in Euclidean relational particle mechanics, as opposed to summing to a constant I in similarity relational 
particle mechanics. In the solution exhibited, the particles expand away from triple collision while relative angle varying oscillations occur, 
involving almost-isosceles to almost-collinear changes in shape.] 

In SSec 6.8's example, i\ and l 2 oscillate boundedly, so neither of these (or Ii or l 2 is a good clock parameter from the 
point of view of monotonicity. There is again some scope for variation in relative angle including 'sporadic' amplitude 
variations. 



6.11 Normal coordinates for scaled triangleland with multi-harmonic oscillator potential 

The working of Sec 5 holds again (using now x, in place of n and iin but radii are unaffected by rotations and so 
cancel out giving the same rotation and 9, $ to 8n, 'I'n coordinate change as before. For Euclidean relational particle 
mechanics, one can uplift from the preceding parts of Sec 6 by inserting the above extra steps into the triangleland case 
of Sec 3.1. 

The potential is now 

{■K?p? N + K^p 2 N }/8l = A N + BncosGn . (150) 

Then the Jacobi-type action for scaled triangleland with general multi-harmonic oscillator potential in shape-scale vari- 
ables is 



S = 2 / dA 



'{i +i 2 {e 2 



sin 2 e N $ N } 



.BncosGa 



(151) 



The Euler-Lagrange equations and energy integral that follow from this are, respectively and after discovering the con- 
served quantity J and eliminating it (|120[ I121|) with N-subscripts appended. The momenta are (|171[) with N's appended 
in the last two and the Hamiltonian and the energy constraint are (|172[) treated likewise. 

I have then rotated the above two exact solutions for the special case to obtain solutions to the general case, but these 
are too lengthy to include in this paper. 
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7 Conclusion 

7.1 Results Summary 

Relational particle models are useful as regards the long-standing absolute versus relative motion debate, and also due to 
structural similarities with the geometrodynamical formulation of General Relativity, for Problem of Time in Quantum 
Gravity. 1- and 2-d relational particle mechanics are tractable due to the simple nature of their configuration space 
geometries: these are, respectively, § fc and CP fe for 1- and 2-d similarity relational particle mechanics. Additionally, for 
the 3-particle case of the latter ('scalefree triangleland'), CP 1 = § 2 holds, making this case even more tractable. This and 
its Euclidean relational particle mechanics counterpart ('scaled triangleland') furbish this paper's particular examples. 

I consider models with general multiple harmonic oscillator type potentials which, as compared with the earlier study 
|17j . include the new feature of relative angular momentum exchange between the two constituent subsystems. I get there 
by first considering the 'special' subcase (which has no relative angle $ dependence in its potential and involves no relative 
angular momentum exchange) and its 'very special' sub-subcase (which has constant potential in one presentation). I 
then identify scalefree triangleland's very special case's mathematics with that of the linear rigid rotor; the special case is 
then analogous to that with a background electric field aligned with its axis. The Euclidean relational particle mechanics 
special and very special cases' mathematics reduces to some of the mathematics that arises in the Kepler-Coulomb 
problem. Finally, I use a rotation or normal coordinates construction to cast the general (^-dependent, relative angular 
momentum exchanging) case as the special case in the transformed coordinates. This has the mathematics corresponding 
to the analogue background electric field being unaligned with the axis. In each case I use the standard spherical or 
planar mathematics that I have been able to cast the problem into so as to obtain solutions, and then map back to 
provide physical interpretation in terms of various mechanically-significant quantities that are more intuitively associated 
with the original relational problems: (barycentric partial moment of inertia, $) variables, mass-weighted Jacobi inter- 
particle (cluster) coordinates, and by sketches, what the particles themselves are doing. In returning to these various 
levels, the standard spherical, planar and flat space mathematics becomes unusual and nontrivial. 

Highlights of my results are as follows. 

1) The very special multiple harmonic oscillator like potential scalefree triangleland problem is straightforwardly soluble 
on the sphere and retains a manageable form in terms of the underlying mechanical variables. 

2) While the special multiple harmonic oscillator like potential scalefree triangleland problem is also classically exactly 
soluble on the sphere, its form is very complicated, so I just provide more manageable large and small asymptotic solutions. 

3) In the stereographic plane, the small asymptotics solution has the mathematics of the 2-d isotropic harmonic oscillator 
[(k, constant) x (radius) 2 potential, including the upside-down (k < 0) and degenerate (k = 0) cases], whereby I obtain 
complete control of it and then characterize it in terms of the underlying mechanical variables. 

4) The large asymptotics solution maps to the small asymptotics solution again, under inversion of the radius, so I also 
obtain complete control of it and can characterize it in terms of the underlying mechanical variables. Moreover, this is 
important beyond the case with harmonic oscillator like potentials, as scalefree triangleland exhibits universal large-scale 
behaviour, various cases of which I can now understand through their being related by the inversion to the various cases 
of conic sections that occur for k x (radius) 2 potentials. 

5) The very special and special multiple harmonic oscillator Euclidean relational particle mechanics problems are also 
classically exactly soluble in the conformally-related flat space in which they repectively take the forms of the Kepler- 
Coulomb problem and a nonstandard composition of parabolic coordinate subproblems thereof. 

6) Then by a rotation or normal coordinates construction, the general relative-angle dependent similarity and Euclidean 
relational particle mechanics multiple HO (type) potential problems are exactly soluble because they map to their special 
cases in the new coordinates. However, these solutions are very complicated in terms of the underlying mechanical 
variables, and as such I mostly only provide sketches of some aspects of their behaviour. 

7.2 Further tractable cases: 4-stop and N-stop metrolands 

Because for scalefree 4-stop metroland (relational particle mechanics of 4 particles in 1-d) the Rj are different and 
related to 0, $ in a different way, physically interesting potentials in this case generally map to different functions of 
these coordinates for the 4-stop metroland interpretation and for the triangleland interpretation. Thus one needs new 
calculations for 4-stop metroland rather than straightforward deduction from this paper's triangleland results. E.g. the 
harmonic oscillator type potential maps to 



for A = {K x + K 2 + 2A' 3 }/8, B = {-K x - K 2 + 2A 3 }/8, C = {K 1 - K 2 }/A. Thus this model admits direct analogues of 
this paper's 'special' (C = 0) and 'very special' (B = C = 0) subcases. 

While, generalizing to scalefree N-stop metroland, the multiple harmonic oscillator type potential maps to 



V = 



A 3 cos ; 



i 2 9 + Aisin 2 6cos 2 <I> + A 2 sin 2 esin 2 $ 
2 



A + Bcos(26) + Csin 2 6cos(2$) 



(152) 




(153) 
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for rip the unit vector in the Euclidean configuration space K n . Scaled N-stop metroland is also of interest [55] . 



7.3 Further work and Applications 

This paper's models remain to be studied from the perspective of dynamical systems [71;. Its principal application at the 
moment is that many aspects of this work carry over to QM in paper II (similarity case) and |37j (Euclidean case), and 
on towards the study of many Problem of Time strategies (in particular emergent semiclassical time and records theory 
[19] I20[ 131] but also conceivably internal time approaches and histories theory, as well as investigation of various Quantum 
Gravity and Quantum Cosmology applications such as the problem of observables, operator ordering, closed-universe 
effects, finite- universe effects and the study of small inhomogeneities/clumps). Some of these applications would benefit 
from study of further models: with other potentials (for which universal large asymptotics results in this paper will be 
useful), the above 4-stop and N-stop metrolands of comparable tractability to this paper, as well as somewhat harder 
relational particle mechanics in 2-d of N > 3 particles (requiring CP N_2 geometry based methods) and 3-d models that 
are much harder [64] even for modest values of N. 



Acknowledgments 



I thank Dr Julian Barbour, Dr Brendan Foster and Miss Anne Franzen for discussions, the two Anonymous Referees 
for comments, Professor Gary Gibbons for references, the organizers of "Space and Time 100 Years after Minkowski" 
Conference at Bad Honnef, Germany, the Perimeter Institute and Queen Mary's Relativity Group for invitations to speak 
and hospitality and Dr Julian Barbour also for hospitality. Peterhouse for funding this work in 2006-08. Professors 
Malcolm MacCallum, Gary Gibbons, Don Page, Reza Tavakol and Jonathan Halliwell, and Dr's Julian Barbour and Fay 
Dowker, for support in the furthering of my career. My Wife for help with assembling the figures, and my Wife, Alicia, 
Amelia, Beth, Emma, Emilie, Emily, Joshua, Luke, Simeon and Will for keeping my spirits up. 

Appendix A: Mechanics on an in general curved configuration space 

For a general finite theory of the quantities Qa with a curved configuration space, consider the Jacobi action 

S = 2 J dAv/T{U + E} (154) 

where U(Qa) is minus the potential energy V(Qa), E is the total energy and T is the kinetic energy, 

T = M AB Q A Q B /2 , (155) 

for A^ab the curved configuration space metric. 



This action works as follows. For * = y{U + E}/T ', the Euler-Lagrange equations are, in geometrical form, 

DM 
dQ~A 



A** +r A BcQ B» Q C* _ _ ( 15 g) 



and there is an energy first integral 

M AB Q A *Q B */2 + V(Q C ) = E . (157) 

Appendix B: Momenta, Hamiltonians and energy constraints 

These are important as regards the passage to Quantum Theory in paper II and [37] . 
B.l General curved configuration space mechanics 

The conjugate momenta are 

Va = M ab Q B \ (158) 
and there is then as a primary constraint the quadratic energy constraint 

H = N ab V a Vb/2 + V= E (159) 

for A/" AB the inverse of A^ab and H the Hamiltonian for the system. This is propagated by (|156p . Such energy constraints 
that have quadratic but not linear dependence on the momenta are analogous to the GR Hamiltonian constraint and 
carries associated with its form the frozen formalism aspect of the Problem of Time [55] [33J . 
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B.2 Scaled triangleland in (ti, l 2 , 3?) coordinates 

The conjugate momenta are 

p; = t* , p $ = + 4} . (160) 

The classical Hamiltonian and energy constraint are then 

h = ^e{to 2 + ^}+ v = e - ( l61 ) 



In the ^-independent case, P$ = J, constant, so one has furthermore 

2 r t2 



" = ig{<^ + f}+V = E 



(162) 



B.3 Scaled triangleland in (Ii, I 2 , «&) coordinates 



These are useful in the context of a ^-independent potential energy in the special and very special cases, for which the 
conjugate momenta are 

Pi = Ij/Ii , P& = 4l!l 2 $7l = J , constant . (163) 
The Hamiltonian and the energy constraint are then 

B.4 N-paricle d-dimensional preshape space theory and scalefree N-stop metroland 

The conjugate momenta are 

P q = |nsin 2 e p -|e| . (165) 
The Hamiltonian and the energy constraint are then 

H = ^£ hiAtt + V(e ^ = E • (166) 

B.5 Scalefree N-a-gonland and the exceptional case of triangleland 

The conjugate momenta are 



Vn * ii + ii^ii 2 {i + ^ii 2 } 2 }^ ' v& * ii + ii^ii 2 {l + m 2 } 2 } 71 * 71 ^' (167) 

The Hamiltonian and the energy constraint are then 

H = 2{1+W|P> {i SPQ + nPnQ } V ^ V ^ + {w p + irj^Pe,} +V(^-,6p) = E . (168) 

For the specific example of scalefree triangleland with harmonic oscillator like potentials in stereographic coordinates, the 
conjugate momenta are 

PK = n* , p@ = n 2 e* (169) 

and the Hamiltonian and the energy constraint are given by 

- If 2 p% \ K{R? + L^cos$ + K 2 _ E 

2\ Pn+ 'R?] + 2{l+^2}3 ~{l+^2}2- [UU > 

While, in spherical coordinates, one has 

p @ = 0* , 39$ = sin 2 0<T (171) 

and 

H = ^ \p% + -^f- 1 + A + PcosG + Csin0cos$ = E . (172) 
2 ^ sin J 
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B.6 Scale-shape formulation of scaled triangleland 



The conjugate momenta are 

pi = I* , p& = I 2 6* , = I 2 sin 2 e$* . (173) 
The Hamiltonian and the quadratic energy constraint are then 

2 2I 2 2I 2 sin 2 8 41 
In the special case, p$ = J . The Hamiltonian and the energy constraint are then 



H = ^ + ^ + ^^ + ^^-E. (174) 



Z_pi 2 ps 2 J 2 v(i,e)_ e 

H _ ^ + 2I 2 + 2I 2 sin 2 e + 41 E ■ (175) 



In particular, for scaled triangleland with multi-harmonic oscillator potential, this is 

2 2 2 p 

^ = pi Pe m A BcqsQ CsiaQcos $ = E = — , (176) 
2 2I 2 2I 2 sin 2 6 41 y J 

which is close to but not exactly the same as the classical Hamiltonian for an atom in a background homogeneous electric 
field pointing in an arbitrary direction. And in normal coordinates, the Hamiltonian and the energy constraint are 

H = 5 + ^% + 9% +A N + B N cose = ^ , (177) 
2 2I 2 2I 2 sin 2 9 41 y J 

which is close to but not exactly the same as the classical Hamiltonian for an atom in a background homogeneous electric 
field pointing in the axial 'z' direction. 

Appendix C: Physical interpretation of this paper's relative angular momentum quantities 

The relative angular momentum quantity in scaled triangleland in simple relational variables is 

J = hh$*/I . (178) 

This is equivalent to the scale-shape spherical polar form 

J = I 2 sin 2 6$* (179) 

by (|134|) and (I109[) . It has the following interpretation. 

JI = I X I 2 $* = hhie; - 8{} = I X L 2 - I 2 Li = {Ij + I 2 }L 2 = -{Ij + I 2 }Li (180) 

where the fourth equality uses the zero angular momentum constraint, and so, as Ii + 1 2 = I, 

J = L 2 = -L x = {L 2 - Li}/2 . (181) 

So it is interpretable as the angular momentum of one of the two constituent subsystems, minus the angular momentum of 
the other or half of the difference between the two subsystems' angular momenta, which is a relative angular momentum 
presentation. 

Using spherical coordinates, scalefree triangleland's relative angular momentum quantity is 

J = sinW = = lWe$ * = I (182) 

I 41 I I V ' 

by 62), (plD and (TTTH]). Thus, by (TTgll . 

J = J/I = L 2 /I = -L x /I = {L 2 - Lx}/21 , (183) 

but I is constant in similarity relational particle mechanics, so this is still, up to a constant of proportion, the angular 
momentum of one of the two constituent subsystems, minus the angular momentum of the other or half of the difference 
between the two subsystems' angular momenta. Additionally, it has the dimensions of rate of change of angle, which 
makes sense since on the sphere only angles are meaningful. 

In either case, ^-independence in the potential corresponds to there being no means for angular momentum to be 
exchanged between the subsystem composed of particles 2, 3 and that composed of particle 1. 

For scaled triangleland, one can consider the configuration space vector 

Isin9cos$ \ / 2tit 2 cos$ \ 
IsinGsin* ] [ 2iit 2 sm$ | 
IcosG I = \ I ■ (184) 
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Then from this and its conjugate momentum P, the vector 

1 = 1 x £ (185) 

can be formed, which is conserved in the very special case; J is the axial 'Z' component of this, Jz, while the new 
components are 

2 

3 * = 2 T 2 {i L i + L 2}{^2 ~ t2ii}sin$ + {i\ - ^}tit 2 cos$$*} (186) 
ij + i 2 

and 

2 

J Y = 2 {{>i + t4}{ti4 - t2ii}cos$ + {l\ - tl}iii 2 sin$$*} . (187) 

h + L 2 

Additionally, the configuration space Laplace-Runge-Lenz type vector 

Q = P x J - EI/I = (188) 

2tit 2 cos$ \ / {t*t2 + t!ii}cos$ - iii 2 sin$$* 

2iit2sin$ ] [ {t*t 2 + i 2 ti}sin$ + i 1 /, 2 cos$<I>* 

2 ~~ L l 

(189) 

which is also conserved. However, this only furbishes one further independent conserved quantity due to the interdepen- 
dences 

J ■ Q = and Q 2 = E 2 - 2AJ 4 /E 2 . (190) 




4{{<-l + 4}Wl + i 2 2 } + l 2 !^* 2 } ~ -2—^2 f \ I -±{ll+Ll}{lHl+l2L*2} \ L 2 L* 2 - Ht\ 



The very special multiple harmonic oscillator like potential case of scalefree triangleland also has not just a conserved 
quantity J but a conserved vector J_ of which J is the axial 'Z' component (so that J_ = J/I). ^-independent scalefree 
triangleland can still have yet more conserved quantities, but these are of a rather more complicated nature along the 
lines described in e.g. [771 175] , 
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